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Note to be inserted at p. 187, Part I. 1882. 



CORRECTION TO A PAPER ENTITLED 

" On the Stresses caused in the Interior of the Earth by the Weight of Continents and Mountains" 

Since this paper has left my hands I have discovered an error in the work. The 
error does not affect the physical conclusions, except in one unimportant respect ; had it 
done so I should probably have found it out long ago. 

Throughout the paper the normal stresses P, Q, R require an additional term W.. The 
only function of these stresses used for obtaining physical results is P — R, and it remains 
unchanged when the correction is made. § 10 must however be erased. 

The error takes its origin in § (1), Thomson's solution (1), when reduced to the form 

applicable to the incompressible solid, is the solution of the equations -• -£ -f ^a 2 = -^-~, and 



J*AJ 



two others. The solution required is that of — ~~ + vtjol 2 = 0, and two others. The W 

involved in my solution is not the potential of a true bodily force, but only an "effective 
potential" producing the same strains as those due to the weight of the continents and 
mountains, but causing a different hydrostatic pressure. When therefore p is determined 
from Thomson's solution, that p is really equal to p+ Wi of the problem of the continents. 

Hence equation (3) should be p=— (l+yjT^i, instead of _p==— yW^. The correction to 

(3) must be carried on through the rest of the paper, and obviously it merely adds Wi to 
the stresses P, Q, R, leaving P — R unchanged. 

The error, would have been avoided had I, as suggested on p. 190, worked directly from 
the equations of equilibrium of the elastic incompressible solid, instead of from Thomson's 
solution. 

When the solid is compressible, this method of " effective potential " [see " The Tides 
of a Viscous Spheroid," Phil. Trans, Part I. 1879, pp. 7 — 9] for including all the effects 
of gravitation, is not applicable without certain additional terms in a, /3, 7. Hence § 10 is 
erroneous, inasmuch as the expressions for the strains and stresses are incomplete. The 
correction of § 10 (which is not difficult) would require too much space to be carried out 
in this note. 

G. H. DARWIN. 

Aug. 1, 1882. 
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In this paper I have considered the subject of the solidity and strength of the 
materials of which the earth is formed; from a point of view from which it does not 
seem to have been hitherto discussed. 

The first part of the paper is entirely devoted to a mathematical investigation, based 
upon a well-known paper of Sir William Thomson's. The second part consists of a 
summary and discussion of the preceding work. In this I have tried, as far as possible, 
to avoid mathematics, and I hope that a considerable part of it may prove intelligible 
to the non-mathematical reader. 
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THE MATHEMATICAL INVESTIGATION. 
§ 1. On the state of internal stress of a strained elastic sphere. 

Let there be a homogeneous elastic sphere, for which co— -\v is the modulus of com- 
pressibility (or ineompressibility, as I shall call it) and v the rigidity.* Take the 
centre of the sphere as the origin for a set of rectangular axes x 9 y, z. Let the sphere 
be subjected to no surface stresses, let it be devoid of gravitation, but subject to 
internal force such that the force acting on a unit volume- of the elastic solid is expres- 
sible by a gravitation potential W h a solid spherical harmonic of the i th degree of the 
coordinates x, y, z. 

Let w be the density of the elastic solid, a the radius of the sphere, and r the radius 
vector of any point measured from the centre of the sphere. 

Sir William Thomson has investigated the state of internal strain produced under 
the conditions above described. If a, /3, y be the displacements his solution is as 
follows :— 



a=(E^-F^)^-Q t ^ + ^(W t r'^ ^ 



where 



Ei= 



Ct'OQ (Xh)g 



2(^-l)v{[2(^ + l) 3 + l>-(2^ + l)v} 



i ~~2{2i + r)v{[2(i + iy + l]co--(2i + l)v} 



V (I) 



n 10) 



(2i + l)v{[2(i + iy-t-l]co-(2i+l)v} 



j 



and similar expressions for (3 and y. t 

Now let P, Q, R, S, T, U be the six stresses, across three planes mutually at right 
angles at the point x, y 9 z, estimated as is usual in works on the theory of elasticity. 

Let P, Q, R, be tractions and not pressures, and let p be the hydrostatic pressure at 
the point x 9 y, z. 

Then P+Q+R being an invariant of the stress quadric, we have, 

if 8=:-^+—- +- ~-_ so that S is the dilatation, then according to the usual formulas,! 

dx ' dy dz 9 & + 

* The phraseology adopted by Thomson and Tait (first edition) and others seems a little unfortunate. 
One might be inclined to suppose that compressibility and rigidity were things of the same nature ; but 
rigidity and the reciprocal of compressibility are of the same kind. If one may give exact meanings to old 
words of somewhat general meanings, then one may pair together compressibility and "pliancy," and call 
the moduli for the two sorts of elasticity the " incompressibility " and rigidity. 

f Thomson and Tait's 'Nat. Phil.,' § 834, (8) and (9) ; or Phil. Trans., 1863, p. 573. 

J Thomson and Tait's « Nat. Phil./ § 698. 
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da 



do: 



= {G) — V)0-f-Zv 

rp /dot dy\ 

jizzzyi — -4- — } 

\dz dxj 

and the other four stresses are expressed from these by cyclic changes of P, Q, R; 
, 1, U ; a, p, y ; x, y % %. 
The first task is to find p. 
Now by adding P, Q, II together we have, 

p=~(co-iv)8 
We must now find 8 from (1). 
By differentiation 

da d?W- dW> d% d 

and similar expressions for d/3/dy, dyjdz. 

Now Wi, WiT~ %i ~~ l are spherical harmonics of degrees i, — i— 1, and are also homo- 
geneous functions of the same degrees. If therefore we add the three expressions 
together, and note the properties of harmonics and of homogeneous functions, we have 

S= — 2iF i W i +(2i+3)(i+l)G i W i 

Omitting for brevity that part of the divisors in the expressions for F and G which 
is common to both, 

— 2iFi+(2i+3)(i+l)Gi=—i(i+l)(2i+fy^ 

= i(2i+l)v 

and we have, on introducing the omitted denominator, 

# 

. § = : 1 W. 

[2(i + lf + l]u-(2i+l)v * 

r~ [2(i+ 1) 3 + 1> —(2*4- l)w * 

Throughout the rest of this paper (excepting in § 10) the elastic sphere will be 
treated as incompressible, so that <w is to be considered as infinitely large compared 

with v. 

Henceforth I write 

JL «— • £J I V \ J- / I A » » » • » « « ■ » » \ iLA ] 

and when o is infinite compared with v, we have, 

* 

if m— i mmmm ff*£ * * * • * * # * > * * \^/ 

Also we may put 
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P 



T: 



, ^ da 
-p+2v~- 



dx 



v 



da d<y 



> 



dz dm. 



i * * » 



• • (4) 



And on putting co Infinite in (I) we have 



a\ 



1 
Tv 



% 



i V + 2) , _ (i + 1) (2i + 3) 1 rf Wj _ 
LL2(*'-1) 2(2*4-1) J die 2i + l' 



„2i+S_ 



d 



-_ "V 



rfe 



(W^r-^) 



and symmetrical expressions for /3 and y. 



• • ( 5 ) 



The hydrostatic pressure might have been found from this general solution for the 
case of incompressibility, but in order to do so it would have been necessary to go back 
to the equations of equilibrium of the solid, and I prefer to deduce it from Sir William 
Thomson's solution in the more general case. 

Since 



-* „ j*_p *|~, -um 

A-+sfL( Wp-*- 1 ) = - (2i+ l)x Wi+r*-p 



dx 



and since 



dx 



(i+'l)(2i+3>) + 2i=(2i+l)(i+S) 



we may write a as follows :- 



a= 



2Iv 






(i+3)?* 3 1 



}dW i 
ax 



-\-2ixWi 



»*»•»• 



(6) 



^^^ $yfi£ £kfl£ fl^f 

In order to find the stresses P, Q, &c, we must now evaluate —> 7 ~> -~> &c. 

fe flfe a#? 

Differentiating (6) with regard to x, — 



da \i(i + 2) 2 

<x^ ^ — 1 



({+3)r 



2 



-^-ea^y-M^W^- (7) 

ax* dx v 7 



Differentiating with regard to z 9 - 



^ -r da 

2lv~ — 

dz 



i(i+2) o /• , ox o 



dxdz 



-\-2] ix 



. dW, 



dz 



(i+S)z 



dWi 
dx 



• a • 



and by symmetry 



9 f *? — 



a 



—. — ~ or— u+3)7^ 

^ — 1 v 7 



cfedfe I dx 



(i-\-3)x 



dWi 

dz 



(8) 



• « » 



(9) 



Adding (8) and (9) together and dividing by 2, we have 



7 /<&* , *y\ 



i(i+2) 



a 2 — (i+3)r 



,2 






— 3 



£C— — + Z- 



d" 



<i# 



(10) 



Hence from (3) (4) (7) and (10) we have. 
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/P= 



IT= 



- — zr-a? 

% — \ 
i{i + 2) 



— ('i*+3)r 2 



<PW; 



a 



(i+Z) 



T* > 



dxdz 



cl Wi 

dx ' 

0/ dWi , dm 



> 



> 



dz 



di 



(14) 



where /=2(i+l)*+L 

The expressions for Q, R, S, U may be written down from these by means of cyclic 
changes of the symbols. 

These are the required expressions for the stresses at any point in the interior of 
the sphere. 

In order to find the magnitude and direction of the principal stress-axes at any 
point it would be necessary to solve a cubic equation. The solution of this equation 
appears to be difficult, but the special case in which it reduces to a quadratic equation 
will fortunately give adequate results. It may be seen from considerations of symmetry 
that if Wi be a zonal harmonic, two of the principal stress-axes lie in a meridional 
plane and the third is perpendicular thereto. Moreover the greatest and least stress- 
axes are those which lie in that plane, and the mean stress-axis is that which is 
perpendicular thereto. If this is not obvious to the reader at present, it will become 
so later. 

1 shall therefore take Wi to be a zonal harmonic, and as the future developments 
will be by means of series (which though finite will be long for the higher orders of 
harmonics) I shall attend more especially to the equatorial regions of the sphere. 



§ 2. The determination of the stresses ivhen the disturbing potential is an 

even zonal harmonic. 

If be colatitude the expression for a zonal surface harmonic or Legendre's 
function of order i is 



cos* 0— ) - A cos*"" 3 sm 2 0+ K -tt—" — cos 2- " 4 sm 4 0— . . . 



4.(1 !) 3 



4 2 (2!) 2 



or if we begin by the other end of the series, and take i as an even number, the 
expression is 



7, 1 



i'3 



sin* 6—jr % sin*"" 2 cos 2 0+ 

2! 4! 



iHi-2f . . 



sin*"" 4 cos 4 



... 



This latter is the appropriate form when we wish to consider especially the equatorial 
regions, because cos is small for that part of the sphere. 

There is of course a similar formula when i is odd, but of this I shall make no use. 
Now let p % z=.y^-\-z % , so that sin 0=p/r, cos 0=z/r. 
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Then we may put 



VV « — P o\P Z ^ A\ P * 



I* 

21 



4! 



6! 



p'- b z b + . . . . 



(12) 



W 7 * is a solid zonal harmonic of degree i\ but r^Wi requires multiplication by a 
factor (—)^!/2 i {^'!} 2 in order to make it a Legen dee's function. 

The factors by which Wi must be deemed to be multiplied in order that it may be 
a potential, will be dropped for the present, to be inserted later. Or we may/ if we 
like, suppose that the units of length or of time are so chosen as to make the factor 
equal to unity. 

Now let 



Po — x > Pz — 2!' "-f" 41 ' ^ 6 — R! 

Then, dropping the suffix to W for brevity, we may write 

w=$ oP i -^p i -w+p i p i -w-p 6 p i -h*+ . . . 



(13) 



• ft 



• (14) 



I shall now find P, Q, P, T at any point in the meridional plane which is determined 
byy=0. 

In evaluating the first differential coefficients of IF we must not put y=0, in as far 
as these coefficients are a first step towards the determination of the second differential 
coefficients. But in as far as these first coefficients are directly involved in the 
expressions for P, Q, P, and T, and in the second coefficients in the same expressions, 
we may put y=0, and thus write x in place of p. 



dp dp dp A - - o 9 1 



dx 



Then 



dz 



clW 

dx 

dW 
dy 

dW 

dz 



=2/ [same series 



] 



=*[ 



2fop?-*+iP 4 ft-W-6l3 6 ft-W+ . . :_ 



In differentiating a second time we may treat p as identical with x, because y is to 
be put equal to zero. Thus 



dx 2 

dif 

d?W 

dz 2 



;i(i- lJ&aT 8 — {i-2)(i- 3)/3 2 af~V+ (i~ 4)(i- 5)/3 4 aT 6 z' 



• » « 



> 



= -L2^- 2 +3.4/8 4 £c 8 '-% 2 -5.6/8 6 ^- 6 z 4 + . . . 



(15) 
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d?W 



d*W 

dxdy 



=xz[_— 



= 0, 



2(i-2)fox i ~ i +i(i-4)P 4! x i -W-6(i-6)/3 6 x i - 8 z i + . . .] 



dydz 







Also treating p as identical with x, and putting y=0, 



. (16) 



dW 

dW_ 
dW 

z — ——2^-^+4,^-^-6/3^-^+ . . . 



>..-.. (17) 



dW , dW 

CvW tv# 



y dx ' dy 



xz {(^ -2AK'- 2 -[(**-2)A-4/3 4 ]^-V 

+[(;_4)&-6/? 6 >rV 

0, 



' ' ') 



dW , rfJT. n 



& dy 



. (18) 



These various results have now to be introduced into the expressions (11) for 

X y V£j JLVj ^J JLj U . 

In performing these operations it will be convenient to write J for i(i+2)/(i—l). 

Also r 2 =p 2 +3 2 =# 2 +2 2 > when y=0. 

From these formulas we see that S=0, U=0; which shows that a meridional 
plane is one of the three principal planes, a result already observed from principles of 
symmetry. 

Now 



JPW_ 

do? 



i(»-i)A^+Wt-i)A-(*' 

-[(.-2)(i- 



&W 



2)(i-s)&y-w 

3)&-(t-4)(t-5)j8J«'-V+ • • • 

(i-2)&-(i-4)(B i y-w+ . . . 



^ = i/3 xt+[i/3 -(i-2)&-]xt-W- 
r ^=_1.2 i 8^-[1.2i8 a -3.4i8 4 ]arV+[3.4)8 4( -5.6i8 fl ]a*-V- . . . 



V 



(19) 



J 



2x~+iW=-i/3 x i +(i-4)fox>-W-{i-8)/3 i x i -W+ . . . 



dx 



?— 4/y4. 



2y~+ i W= ifaf - * A^" V + i&rt 



» • • 



> . . (20) 



-2^f +i WsnA/rf- (i- 4)/3 2 ^-% 3 + (i- 8)^-%*- 



cfe 



*s 



MDCCCLXXXIL 



2 c 
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r ^ z =^{-^-^)/3 2 o ( f-"-[2(i-2)/3 2 -i(i-4:)$ i y-^ 

+[4(*-4)/3 4 -6(;-6) y 8 6 >'-% 6 - . . .} . . (21) 

Then multiplying (19) by — ({+3), (20) by 3, and (15) by Ja?, and adding them 
each to each, we get the expressions for P, Q, R. 

Also multiplying (21) by —({+3), (18) by —3, and (16) by Ja % and adding, we get 
the expression for T. The results are 

7P=-[(i+3)i(»-l)+3t]/8 af 

+ {{(i+3)(i-2)(i-S) + ^i-A)}^-(i+B)i(i-l)^y-^ 

-[{(i + 3)(i-4)(i-5) + 3(i-S)}^-(i+3)(i-2)(i-3)^']x i -^ 

+ [{i+S)(i-G)(i-7) + S(i-12)}^ 6 -(i+3)(i-4:)(i-5)^']x i - 6 z ii - . . . 

+Ja 2 [^i-l)^- 2 -(^-2)(i-3)^-%H(^-4)(*'-5)^- 6 2 4 - . . .] 

/R=[( i >3)a.2 / 82+3^ ]^-[(?;+3).3.4y8 4 -{(i+3).1.2-3^'-4)} / 8 2 > ! '- 2 3 2 
+ [(i'+3).5.6^ 6 -{({+3).3.4-3(i-8)}jS 4 ,> 8 '-% 4 
-[(i+3).7.S/3 s — {(i+S).5.6 — S(i— 12)/8 6 >'- R z 6 + . . . 
— J"a 2 [1.2^'- 2 -3.4^'-% 2 +5.6 / 8 6 a;»- 6 z 4 — . . .] 

7Q=-[(»43)i-3i]8 !C i 4.[|(i+3)(i-2)-3ij J e 3 -(*'+3)i i 8 ]af'-V 
-[{(»+3)(t-4)-3i}/8 4 -(i+3)(i-2)i8 a ]a!'-V 
+ [{(i+S)(i-6)-3i}^-(i+3)(i-4 : )/3 i ]x i -V- 

+ JaZ[ix i -' 2 -(i-2)fax i - i z i +(i—i)p 4! x i - (i z i — . . 

f=[{(^+3)2(^-2)+3.2}A-3^ > 4 - 3 

-[{(^3)4(t-4) + 3^ 

+[((H3)6(i~6) + 3,6}ft»{(i+3)4(^4) + 3(i--4)}^r¥ 
^[{(^3)8(^8) + 3.8}/3 8 ~{(^3)6(^6) + 3(i«6)}^ 6 ]^'-^ 6 + . . . 
-Ja 2 [2(i~2)/3^ . . .] 

The general law of formation of the successive coefficients is obvious, and it is easy 
to write down the general term in each of the eight series involved in these four 
expressions ; the best way indeed of obtaining the formulas given below is to write 
down and transform the general term. 

The semi-polar coordinates used hitherto are not so convenient as true polar 
coordinates ; I therefore substitute r, radius vector, and I, latitude, for the x, % 
system, and putting x=r cos 1, z=r sin I write 



• • • 
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P=r { cos i l(A +A< i timH+A 4l t&XL 4, l-\- . . .) 

+aV- 3 cos i-H[B +B z tan H+B^ tan H+ . . .) 
K=r' cos f?( C + C a tan 3 Z+ C 4 tan H+ . . . ) 

+aV~ 3 cos *- 3 Z(Z) +Z> 3 tan 3 Z+D 4 tan 4 Z+ . . .) . 
T =r» sin Z cos - 1 Z(.E + J £> tan H+E± tan *Z+ . . .) ■ •* ^ 

+aV-» sin Z cos ^(i^+i^j tan 3 Z+/^ tan 4 Z+ . . .) 
Q = H cos <Z(£ + (t 2 tan 3 Z + # 4 tan 4 Z + . . .) 

+aV" s cos 8 '- 3 Z(iT + JI 8 tan S Z+H 4 tan H+ . . .) 

Then introducing for J" and for the fi's their values in terms of i, I find that the 
coefficients A, B, &c, are reducible to the forms given in the following equations : — 

7^ ==-i{ 4 -(t+2)(t-0)-3.0.(i+l)}+0.(i+3)(t+2)(t+l)=-» s (t+2) 



IA 3 
IA 



0! 
^{((t_0)(t-2)-8.2(.-l)}-i(»+3)(*+0)(t-l) 



2! 



0! 



_ t»(i-2) 



4 — 



4! 



{»( l -_2)(i-4)-3.4(i-3)} + ^(i'+3)(t-2)(i-3) 



2! 



^_ g(j ggt 6/ { . ( ._ 4)( ._ 6) _ 36( ._ 5)} _z ! (j-^ ( . + 3)( ,_ 4)( ,_ 5) 



&c.=&c. 



I&= 



0' 



iA=- 



'2 






Zfl. 



4" 



&C.= 



~i^ 4! L (»- 4 )(— 6 ) 

&c. 



(23) 



/C = i{(»+8)^-P(l(-2)-l)+8.6.l]}=t[(»+l)(»+2) + l] 






0! 



-/Co — 



«2 



2 



1L/A, — 



4 



— 2i{(*+3)(t— 2) 3 -[>*(3.0 — 1)+3.2.3]} 

^~ 2) %+3)(Z-4) 3 -|>-(5.2-l) + 3.4.5]} 



ia = 



#(i-2)\i-&Y 



6 



6! 



{(i+3)(**-6) 3 -[i(7.4-l)+3.6.7]} 



&c.=&c. 



JZ>n = 



t'(* + 2) i 3 



0' 



ID„ = 



2' 



Zfr = - 



4, 



&C.= 



i-1 0! 
i(i+2) i*(i-2y 
i—1 2! 
*(* + 2) t^(t-2) g (t -4)» 
T-T 4! 

&c. 



(24) 
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Zfin = 



0' 



jyC 1 ( 3 *+ 3 )-*(*-°)(*+ 1 )]=i(»+ 1 )(* 8 8 ) ^ 



J^=^(.-2)[8(5t+3)-t(t-2)(t-l)] 



3! 



IE 



4' 



^'- 2 ) 2 (,-_4)[5^ + 3)_^_4)(,-_ 3 )] 



6! 



!&= 



6 



^ _2)2( '-%-6)[7(9i+3)-^-6)(*~5)] 



7! 



&e.=&c. 



if: 



0' 



1F = 



t *(t»-4) 
■i — 1 
^_ 4 ) (^_2)(i— 4) 



2' 



IF,.= - 



i-1 3! 

#(^2_4) (^-2)(i:-4) 3 (i-6) 



4 



/F»= 



i-l 5! 



6 



t-1 

&c.=&c. 



7! 



r 



-> 



• (25) 



1 



0! 



IGc 



;-z 



^ {i{ i-2)-3.2}-^(i+3)i 



2! 



0! 



10. 



4 



^*- 2 ^ t -(t_4)_8.4}+l(t +8 )(*. 



4! 



2 7 



2) 



&C.=&C 8 



IHr\ = 



0" 



if* + 2) 1. 
t-l 0! 



Jl JlJ-c 



2' 



in 



^ + 2) ^-2)* 
i-1 4! * ' 



&c.=&c. 



• (26) 



These sets of coefficients are all written down in sucli a form that the laws of their 
formation are obvious, and the general terms may easily be found. I have computed 
their values from these formulas for the even zonal harmonics of orders 2, 4 ; 6, 8, 10, 
1 2 ; the results are given in the following tables both in the form of fractions and of 
decimals approximately equal to those fractions. 

The G*s and H'b were not computed because their values were not required for 
subsequent operations* 
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Table I. — The coefficients for expressing the stress P. 



• 


A Q 


A 


■^4 


A 6 


So 


B, 


^ 


^ 


2 
4 


16 

1 9 

-•8421 


22 

19 

-1-1579 






_1_ l 6 

'19 

+ -8421 








32 

1 7 

-1-8824 


4-2.8 
"51 

+ -5490 


J.48 

+ 2-8235 




4_ 3 2 
'17 

+ 1-8824 


. 12 8 

5 1 

-2-5098 






6 


32 

1 1 

-2-9091 


+ 18 
+ 18-0000 


_1_ 1 84 

"""11" 

+ 16-7273 


27 2 

~5 6~ 

-4-9455 


4-32 
'11 

+ 2-9091 


JLl 52 

5 5 

-20-9455 


J. 2 56 
l "55 

+4*6545 




8 
10 


6 4-0 

16 3 

-3-9264 


_1_ 1 O 3 2 S 
' 16 3 

+ 63-3620 


2 112 

"16 3" 

-12-9571 


12 16 

16 3 

- 74-6012 


4-6.4.0 

'163 

+ 3-9264 


— 7 6800 
114 1 

-67-3094 


1 9 2 1 6 
' 114 1 

+ 80-7713 


8 19 2 

114 1 

-7*1797 


4 00 

8 1 * 

-4*9383 


_L 3 6 1 3 O 
' 2 4 3 

+ 148-6831 

4-12316 
' 113 

+ 285-9823 


6 9 2 00 

2 43 

-284-7737 


5 6 000 

2 4 3 

-230*4527 


4-400 

' 8 1" 

+ 4-9383 


112000 

"7 29 

-153-6352 


13 2 OH 
' 7 2 9 

+ 438-9561 


5 12 0. 

2 4 3 

-210-6996 


12 


6 7 2 

113 

-5-9469 


13 8000 


12 4 


4-67 2 

'113 

+ 5-9469 


3 6 2 8 8 


4_J,8S 1600. 
' "124 3" 

+ 1513-7570 


.g.1 50400 

"12 43 

-1730-0080 


113 

-1221-2389 


"'113 

+ 212-3894 


12 43 

-291-9389 



Table II. — The coefficients for expressing the stress R. 



4 ■ 

% 

2 
4 


Oo 


o, 


p 


o 6 


D 


^ 


D* 


D» 


4-2.6. 

'19 

+ 1-3684 


+ 1-6842 






3 2 

19 

-1-6842 








_1_ 1 24 

+ 2-4314 


112 

"5 1 

-2-1961 


2 5 6 

"5 1 

-5-0196 




12 8. 

5 1 

-2-5098 


+ 5-0196 






6 
8 

10 
12 


4-3.8. 
~11 

+ 3-4545 


-24 
-24-0000 


— _ 2 8 
"l 1 

-18-9091 


15 12 

1 5 5"" 

+ 9-3091 


J 92 

5 5 

-3-4909 


4.1536 
l 5 5 

+ 27-9273 


5 12 

5 5 

-9-3091 




_1_ 7 2 8 
'16 3 

+ 4-4663 


12 3 5 2 

1 6 3' 

-75-7791 


±4224 

T-"T-8"3" 

+ 25-9141 


_1_ 7 6 2 8 8 

~ 8 15 

+ 93-6049 

4-ioiM 

~ 2 41 

+ 247-5720 


— 5120 
114 1 

-4-4873 


19 2 16 
' 114 1 

+ 80-7713 


2.2 2 88J2. 

114 1 

-107-6950 


J_l 6384 
' 114 1 

+ 14-3593 


113 3 
'2 43" 

+ 5-4733 


4 12 

2 4 3 

-169-5473 


1 S 4 8 
' 24 3 

+ 348*9712 


«— .4 000 
7 29 

-5-4870 


4-JL_2_8_00_Jl 
' 7 29 

+ 175-5830 


J 84000 

"7 2 9 

-526-7490 


4_.20480.p_ 
' 7 29 

+ 280-9328 


17 3 2 
'113 

+ 6-4779 


3 5 8 5 6 
113 

-317-3097 


_1_ 1 5 8 4 00 
~ 113 

+1401-7700 


18400 

113 

-339-8230 


80 64. 

1243 

-6-4875 


4- i0 3 2 0,tt 
» 12 4 3 

+324-3765 


Jg 15 0400 

1243 

-1730-0080 


4.1_5 8048 
' 1243 

+ 2076-0097 
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Table III. — The coefficients for expressing the stress T. 



# 


U 


jE/ 2 


-^4 


JLJ £* 


^0 




F 4 


- 1 6 


2 


'19 

+ •3158 
















4 
6 


_1_ 2 60 
• "5 1~ 

+ 5-0980 


i.80 
•17 

+ 4-7059 






2 6 6 

5 1 

- 5-0196 








+ L4i 
+ 14-0000 


r r 


10 8 

*~*5 5 

-18-3273 




7 68 

~5 5 

-13-9636 


110 2 4 
• 5 5 

+ 18-6182 






1 1 
-5-0909 


8 
10 
12 


4- _* 3 9 2 

T 163 

+ 26-9448 


13 248 

16 3 

-81-2761 


10 3 6 8 

""163 

-63-6074 


.12 44224 
i 5 7 5' 

+ 42-8088 


3 7 2 

114 1 

-26*9238 


J_ 122880 
' 1141 

+ 107-6950 


4 9 15 2 

114 1 

-43-0780 




4_ 1 6 7 O 

' 243™ 

+ 43-9095 


.„_ 7 4 sno 


_i_ 1 7 6 00 


15 7 7 2 8 
' 170 1 

+ 339-3769 


..._ 3 2 


_L. 2 5 600.0 
• 7 2 9 

+ 351-1660 


1 O 2 40 n 


4_ 409 6 
* 5 103 

+ 80-2664 


2 43 

-307-8189 


» ~**2 43*"' 

+ 72-4280 


7 2 9 

-43-8958 

80 640 


2 4 3 

-421-3992 


J.7332 
' 113 

+ 64-8850 


, .0 04 80 


113 7 280 
' 113 

+ 1214-8673 


19 9 840 
' ""113 

+ 883-5398 


_|_ 1075200 
' 12 4 3 

+ 865-0040 


2 5 804 80 

12 4 3" 

-2076-0097 


L1474560 
• 124 3 

+ 1186-2912 


113 

-800*7089 


12 43 

-64-8753 



If W be a 2nd, 4th, or 6th harmonic these tables give the complete expressions for 
P, R, and T ; if W be an 8th harmonic the only further coefficients required are A s 
and Og. 

For the cases of the 10th and 12th harmonics the values in the tables are sufficient 
to give the stresses approximately over a wide equatorial belt, because the series for 
P, R, T proceed by powers of the tangent of the latitude, and the omitted terms 
involve high powers of that tangent. It would hardly be safe however to apply the 
formula — at least as regards the 12th harmonic — for latitudes greater than 15°, 
because the coefficients are large. 



§ 3. On the direction and magnitude of the principal stresses in a strained elastic solid. 

Let P, Q, R, S, T, U specify the stresses in a homogeneously stressed and strained 
elastic solid. Let I, m, n be the direction cosines of a principal stress axis. 

The consideration, that at the extremity of a principal axis the normal to the stress 
quadric is coincident with the radius vector, gives the equations 

(P-\)l+TJm+Tn=0 
TJl+(Q-\)m+Sn=0 
TZ+Sm+(R— X)rc=0 

These equations lead to the discriminating cubic for the determination of X, and the 
solution for Z, m, n is then 

1? m 2 n 2 



(Q-X)(R-X)-S 3 (P~\)(R-\)-T 3 *~(P-X)(Q-~X)-IJ< 
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In the case considered in the preceding sections S and U vanish, and the cubic 

reduces to the quadratic 

(P-X)(R-X)-T 3 =0 
of which the solution is 

2\=P+B±v / (P--R) 3 +4T 3 

m is obviously zero and I, n are determinable from 

P(P-X)=n*(R-X) 

Let 

1= cos 3-, n= sin 3- 

Then it is easily proved that 

cot 25=?-~? . (27) 

This equation gives the directions of the principal stress-axes. 
The two principal stresses N 1? N 3 are the two values of X, so that 

N 1 =i(P+R)+ix/(P-R)^+4Tn 

N,=i(P+B)-iv/(P-R)»+4T»j 

and the third principal stress, which we suppose intermediate in value between N x 
and N 3 , is of course Q. 

When an elastic solid is in a state of stress it is supposed, in all probability with 
justice, that the tendency of the solid to rupture at any point is to be estimated by 
the form of the stress quadric. At any rate the hypothesis is here adopted that the 
tendency to break is to be estimated by the difference between the greatest and least 
principal stresses. For the sake of brevity I shall refer to the difference between the 
greatest and least principal stresses as " the stress-difference," This quantity I shall 
find it convenient to indicate by A . 

We may also look at the subject from another point of view:— It is a well-known 
theorem in the theory of elastic solids that the greatest shearing stress at any point is 
equal to a half of the stress-difference. It is difficult to conceive any mode in which 
an elastic solid can rupture except by shearing, and hence it appears that the greatest 
shearing stress is a proper measure of the tendency to break. This measure of ten- 
dency to break is exactly one-half of the stress-difference, and it is therefore a matter 
of indifference whether we take greatest shearing stress or stress-difference. For the 
sake of comparison with experimental results as to the stresses under which wires and 
rods of various materials will break and crush, I have found it more convenient to use 
stress-difference throughout ; but the results may all be reduced to shearing stresses 
by merely halving the numbers given. 

From (28) we have then 

A = \/(P-R) 8 +4T» (29) 

and the greatest shearing stress at the same point is ^ A ♦ 
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§ 4. The application of previous analysis to the determination of the stresses 

produced by the weight of superficial inequalities. 

I have in a previous paper shown how Sir William Thomson's solution for the 
state of internal strain of an elastic sphere subject to bodily forces, but not acted on 
by any surface forces, is to be adapted to the case of a spheroid (whose small 
inequalities are expressed as surface harmonics) of homogeneous elastic matter, endued 
with the power of mutual gravitation. # Thomson's solution is of course directly 
applicable for finding the state of strain due to a true external force, such as the tide- 
generating influence of the moon, but this forms only a part of the complete solution 
when the sphere has the power of gravitation. He introduced the effects of gravita- 
tion synthetically, but for my own purposes I prefer the analytical method pursued in 
my paper above referred to. 

Suppose that r=a+a"i be the equation to an harmonic spheroid of the i th order, 
forming inequalities on the surface of the sphere, whose density is w. 

Then the causes producing a state of stress and strain in the mean sphere of radius 
a are, first a normal traction per unit area of the surface of the sphere equal to —gwo-i, 
when g is the value of gravity, and secondly the attraction of the inequalities <x^ acting 
throughout the whole sphere. 

The first of these causes (viz. : the weight of the mountains or continents) is shown 
in my paper to produce the same state of strain as would be produced in the sphere, 
now free from surface action, by bodily forces corresponding with a potential 
—gw(r/aycri. 

As regards the second of these causes (viz. : the attraction of the mountains or 
continents), the potential of the layer of matter cri on any internal point, estimated per 
unit volume, is 3gw(r/a) i a'i/(2i+l)A. 

Then adding these two potentials together, we see that the surface inequalities <x ? ; 
produce the same state of strain as would be caused by the bodily forces due to a 
potential — 2(i—I )giv(r/a) i cri/(2i-\'l) > and the surface of the sphere is now subjected to 
no forces. 

This expression is a solid harmonic of the i th degree, and therefore the analytical 

* " On the Bodily Tides of Viscous and Semi-elastic Spheroids, &c," Phil. Trans., Part I., 1879, p. 1 
(see § 2). This paper treats of the state of flow of a viscous sphere, but the problem is exactly the same 
as that concerning elasticity here considered. It is easy to see that if a viscous sphere be deformed into 
the shape of a zonal harmonic, the flow of the fluid must be meridional, and from this we may conclude 
that in the elastic sphere the plane of greatest and least principal stresses is also meridional. This has 
been already assumed to be the case in the present paper. 

f If we could suppose a sphere to have homogeneous elasticity but heterogeneous density, this manner 
of building up the effective disturbing potential would have to be somewhat modified. Such an hypo- 
thesis is somewhat absurd, and I shall regard the sphere as homogeneous. In application to the case of 
the earth I shall however pay attention to the smaller density of the superficial layers by halving the 
height of the actual continents and the depth of the actual seas. 
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results of the preceding sections are directly available for finding the state of stress 

due to continents and mountain chains. 

We must in fact put 

2(*-l) fr\i 

Now in the previous developments the factors involving g, w, &c., have been 
omitted and. Wi has been put equal to a zonal harmonic which had the value unity at 
the equator. 

If we write 

Si=sin*0— — sin'-" s # cos 2 #+&c (30) 



where 6 is the colatitude, and. put h as the height above the mean sphere of the 
elevation at the equator, then hsi=cri and 

Wi was in (12) put equal to r% 

Thus in order to apply the preceding results to finding the stresses caused in a 
sphere, possessing the power of gravitation, by the weight and attraction of surface 
inequalities expressed by 

/ — Cl ~p tl9i » * » . . . . ♦ * • « (da) 

we must multiply the preceding results for P, R, T, Q by 

__2(i-l) gwh , . 

2i + l of •••*•••••• * K 6 *) 



§ 5. The state of stress due to ellipticity of figure or to tide-generating forces. 

When the effective disturbing potential Wi is a solid harmonic of the second degree 
the solution found above will give the stresses caused by oblateness or prolateness of 
the spheroid. It will of course also serve for the case of a rotating spheroid with 
more or less oblateness than is appropriate for the equilibrium figure. 

When an elastic sphere is under the action of tide-generating forces the disturbing 
potential is a solid harmonic of the second degree, and therefore the present solution 
will apply to this case also. 

The formula for the stress-difference admits of reduction to a simple form when 
■i=2. 

On substituting colatitude 6 for latitude I, (22) gives 

P_R=r» sin *ff[(A Q -C )+(Az-Ct) cot *6] + a\B - D Q ) 
T=r* sin# cos 6(E ) 

MDCCCLXXXIL 2 D 
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Then substituting for A, B, C, D, E their values from Tables L, IT., III., and 
making some simple reductions, we find 

P-R=- 1 %[8(a s -r 3 )-r s cos 20] 

2T=T^^sm2^ ** ^ 

Therefore by (29) in the present case, the stress-difference 

A =-& ^Wc^r^f + r^T67^(¥^^y^2e (35) 

In order to find the actual value of A in any special case we shall have to multiply 
(35) by appropriate factors ; the factors will be determined below. For the present it 
will be more convenient to omit the factor -fy, and to reintroduce it along with these 
other factors. 

The formula (35) enables us to determine the distribution of stress-difference 
throughout the sphere in the cases to which this section applies. 

The curves of equal stress-difference are given by the equation 

64(a 3 — r 3 )+r*— 167* 3 (a s — r 2 ) cos 2#=a constant 

The stress-difference at any point on any equatorial radius (for which d=^Tr) is 
clearly given by 8a 3 — 7r s , and along the polar radius (for which #=0) by 8a 3 — 9r 3 . 

From this result it is clear that the stress-difference vanishes at that point on the 
polar radius for which r=-f&\/2=== "9425a ; this is the only point within the whole 
spheroid for which it vanishes. 

When r=a the stress-difference is equal to a 3 , from which we obtain the remark- 
able result that the stress-difference is constant all over the surface. When r=0, 
it is equal to 8a 3 , which is eight times the surface value. 

By means of arithmetical calculation and graphical interpolation I have drawn 
fig. 1 ; it shows the curves of equal stress-difference throughout a meridional section 
of the spheroid. The numbers written on the curves represent the values of stress- 
difference when the radius of the sphere is unity and when the factors above referred 
to are omitted (see Plate 19, fig. 1). 

The point marked is that in which the stress-difference vanishes. Round this 
point are drawn two dotted curves along which the stress-differences are \ and § 
respectively. The remainder of the curves are drawn for equidistant values of stress- 
difference, and are marked 1, 2, 3 ... 8. The curve marked 1 is singular, for the 
whole of the surface forms one branch of it, whilst there is another branch which runs 
below the surface from the polar axis and then rises to the surface at the point where 
cos 20=— -|> that is to say, in lat. 41° 25'. Near the centre the curves are approxi- 
mately circular, and they become somewhat like ellipses as we recede from the 
centre. 
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If this figure be made to rotate about the polar axis, the several curves will of 
course generate the surfaces of equal stress-difference throughout the sphere. 

Writing 3 for the inclination of one of the principal axes to the equator, we have 
by means of the formula (27) 

cot 23-— -^-^=8 \(^\ — 1 1 cosec 20— cot 20 

It would be easy to trace out the changes of direction of the principal stress-axes 
throughout the sphere, but I will only now make the observation that all over the 
surface they are parallel to and perpendicular to the surface, and that at the centre they 
are polar and equatorial, the stress-quadric being of course an ellipsoid of revolution. 

We have next to find the actual amount of stress-difference which arises from 
given ellipticity of form of the spheroid. Putting {=2 in (30) we have 

9i = sin 3 0—2 cos 2 0=S[}— cos 2 0] 

The equation to the spheroid is 






1 + 3 /? '(i - cos 2 0) 



a 



= a[l+e(g— cos 2 0) 



Thus 3A/a is the ellipticity of the spheroid, which we may put equal to e. 
. Then it was shown in (33) §4 that the results for the stresses P, Q, R, &c., are to 
be multiplied by —%gwh/a* 9 and this will of course be also the factor for the stress- 
difference A. 

Then substituting e for 3h/a, and introducing the factor -®$ , which has been omitted 
in considering the distribution of stress within the spheroid, we see that ellipticity e 
gives a stress-difference represented by 

a = ->. V«<>-©T + ('»)- |6 © ! ( 1 -©> os2 " 

If we estimate the forces in gravitation units the factor g must be omitted. 

The expression under the square-root sign is equal to unity at the surface, and to 
8 at the centre. 

Thus the stress-differences, in gravitation units of force, at the surface and at the 
centre are -£%ewa and %jewa respectively. 

To apply this to the case of the earth, take a=637x 10° cm., and w=5'66, and we 
find the surface and central stress-differences to be respectively 152e and 1214e metric 
tonnes per square centimeter. 

If these numbers be multiplied by f3*34, we get the same result expressed in tons 
per square inch. Thus in British units these two stress-differences are 926 e and 

7698e. 

2 d 2 
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If then the ellipticity e be y oVdth, the surface and central stress-differences will be 
nearly 1 ton and nearly 8 tons to the square inch respectively. 

From the Table VII. in § 9 it will appear that cast bi-ass ruptures with a stress - 
difference of about 8 tons to the square inch. 

Thus a spheroid, made of material as strong as brass, and of the same dimensions 
and density as the earth, would only just support an excess or deficiency of ellipticity 
equal to y^tli, above or below the equilibrium ellipticity adapted for its rotation. 

The following is a second example: — -If the homogeneous earth (with ellipticity jij) 
were to stop rotating, the stress-difference at the centre would be 33 tons per square 
i nch. 

Now suppose the cause of internal stress to be the moon's tide-generating influence, 
and let m= moon's mass, and c= moons distance. 

Then the potential under which the earth is stressed is — -f (m/c 3 )(i!r— cos 2 0)w^ 9 or 
according to the notation of § 4 — ^(m/c 3 )wi\ 2 . 

If we took into account the elastic yielding of the earth and the weight and 
attraction of the tidal protuberance, this potential would have to be diminished. To 
estimate the diminution we must of course know the amount of elastic yielding, but 
as there is no means of approximating thereto, it will be left out of account. 

Then it is obvious that the factor by which A, as given in (35), must be multiplied 

in order to give the stress-difference is \mwlc*. Thus the surface stress-difference is 
6 .i 



19 2 



(m/c*)wa* in absolute force units. 



Then putting M for the earth's mass, and dividing by gravity, we have 
•fy(rna?/Mc*)wa as the surface value of A in gravitation units. The central value of 
A is of course eight times as great. 

With the numerical data used above, tw=3605 metric tonnes per square cm., and 
m/M=-^ 9 a/c= e 1 ^ whence the surface stress-difference is 32 grammes, and the central 
stress-difference 257 grammes per square centimeter. 

But this conclusion is erroneous for the following reason. If we suppose the moon 
to revolve in the terrestrial equator, and imagine that the meridian from which longi- 
tudes are measured is the meridian in which the moon stands at the instant under 
consideration, then the tide-generating potential is — -f (m/c 3 )r 2 [-^-- sin 2 cos 2 <£] ; this 
expression may be written f (m/c 3 )r 2 (-|— cos 2 #)+§(m/c 3 )r 2 sin 3 6 cob 20. The former 
of these terms produces a permanent increase of the earth's ellipticity, and is confused 
and lost in the ellipticity due to terrestrial rotation, and can produce no stress in the 
earth. The second term is the true tide-generating potential, but it is a sectorial 
harmonic, and I have failed to treat such cases. Now the first of these terms causes 
ellipticity in a homogeneous earth equal to (f<V#)(fm/c 3 ) according to the equilibrium 
tide-theory. This ellipticity is equal to *1039X 10~ 6 , an excessively small quantity. 
If however this permanent ellipticity does not exist (and the above investigation in 
reality presumes it not to exist), then there will be a superficial stress-difference equal 
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to 152 X* 10 39X1 0~ 6 metric tonnes per square centimeter, and a central stress- 
difference of eight times as much. 

Since a metric tonne is a million grammes this surface stress-difference is 16 grammes, 
and the central 128 grammes per square centimeter. These stress -differences are 
exactly the halves of those which have been computed above. Thus the remaining 
stress-difference which is due to the moon's tide-generating influence is 16 grammes 
at the surface and 128 grammes at the centre per square centimeter. 

A flaw in this reasoning is that stress-difference is a non-linear function of the 
stresses, and therefore the stress-difference arising from the sum of two sets of bodily 
stresses is not the sum of their separate stress-differences. 

I conceive however that the above conclusion is not likely to be much wrong. 

These stresses are very small compared with those arising from the weights of 
mountains and continents as computed below, nevertheless they are so considerable 
that we can understand the enormous rigidity which Sir William Thomson has 
shown that the earth must possess in order to resist considerable tidal deformations 
of its mass. 

§ 6. On the stresses due to a series of parallel mountain chains. 

Having considered the case of the second harmonic, I now pass to the other extreme 
and suppose the order of harmonics i to be infinitely great, whilst the radius of the 
sphere is also infinitely great. 

The equatorial belt now becomes infinitely wide, and the surface inequalities consist 
of a number of parallel simple harmonic mountains and valleys. 

If i be infinitely large, we have from (12) 



W*= 



}-m+h(5J-*°- 



Now let f be the depth below the surface of the point indicated in the sphere (now 
infinitely large) by x, y, z. 

As the formulas given above apply to the meridional plane for which y=0, we 
have p=a— £ 

"Now let b=a/i, then when both i and a are infinite 



p«=a*(l— ~Y==aV~^ 



and since in the limit p/i=a/i=b 9 

= (a»)€~* /6 cosy 
This expression for W involves the infinite factor a\ and in order to get rid of it we 
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must now consider the factor by which it is to be multiplied, in introducing the height 
of the mountains and gravity. 

This factor is computed in § 4 ; it is there shown that if r=a-\-h$j be a harmonic 
spheroid, the factor is — 2(i—l)gwh/(2i-{-l)d l . 

Now if the harmonic i be of an infinitely high order, Si becomes simply cos z/b, and 
the equation to the surface is 

£=—h cos 7 

£ being measured downwards. Thus the harmonic spheroid hsi now represents a 
series of parallel harmonic mountains and valleys of height and depth h, and wave- 
length 27T&. 

The factor becomes —gwh/af, when i is infinite. 

Thus the effective disturbing potential W, which is competent to produce the same 
state of stress and strain as the weight of the mountains and valleys, is given by 



W= —givhe~Z ,b cos j . . . . . , 

Now revert to the expressions (11) for the stresses. 

When i is infinite 7=2^, and they become, on changing x into (a— £) 



(36) 



P 



V-r^+^^+l:^ 



2% 



T= 



df? ' 2v> rff ' 2% 
1, , „«PJT 3 /. ~dW dW 



2f 



'd£fo 2t 8 \ 



ffe 



rff 



Now as shown above a 3 — r 2 =2a£ and a/i=b in the limit; making these substitu- 
tions, and dropping the terms which become infinitely small when i is infinite, we 
have 

d% W d % W ^ 

and by a similar process 



d%dz 



d$W 

&=&>-£, Q=0 



r • i * * i i « . \*J $ J 



d. 



# 5 



-> 



Then from (36) and (37) we have 

P 



*- 



•gwhfe s,v cos 7 



#* 






R=r gwh-e 



-#& 



b 



cos 



T 



& 



£ 



JT • • • • • 



(38) 



t 

q-ivh-e $/b sin r 



Since the stress-difference 



we have 



_> 



ts 



A = V / (P-R) S +4T 



&=2gwhje m . . , , , f , , , . . 



(39) 
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The directions of the stress-axes are given by 



P — K z 

cot 2&=-r=r"= c °t j 



so that 



< ~"~* 27 * • • * * * * * * • • V / 



* 



Equation (39) gives the stress-difference at a depth £ below the mean surface, and 
is very remarkable as showing that the stress-difference depends on depth below the 
mean horizontal surface and not at all on the position of the point considered with 
reference to the ridges and furrows. 

Equation (40) shows that if we travel along uniformly horizontally through the 
solid perpendicular to the ridges, the stress-axes revolve with a uniform angular 
velocity. 

They are vertical and horizontal when we are under a ridge, and they have turned 
through a right angle and are again vertical and horizontal by the time we have 
arrived under a furrow. 

Since the function xC x is a maximum when x=l, the stress-difference A is a 
maximum when g=b,— that is to say, at a depth equal to 1/2tt of the wave-length — 
and is then equal to 2gwhe~ l or in gravitation units of force to '736 wh. It is inte- 
resting to notice that the value of this maximum depends only on the height and 
density of the mountains, and does not involve the distance from crest to crest. The 
depth at which this maximum is reached depends of course on the wave-length. 

Plate 19, fig. 2, shows the distribution of stress-difference, the vertical ordinates 
represent stress-difference, and the horizontal ones depth below the surface. The 
numbers written on the horizontal axis are multiples of h; the distance OL on this 
scale is equal to 6*28, and is therefore equal to the wave-length from crest to crest, 
and the distance OH is the semi-wave-length from crest to furrow. 

In the case of terrestrial mountains w is about 2*8, and if we suppose h to be 2000 
meters, or a little over 6000 feet, we have the case of a series of lofty mountain chains 
—for it must be remembered that the valleys run down to 2000 meters below the 
mean surface, so that the mountains are some 13,000 feet above the valley-bottoms. 

Then h=2 X 10 3 , w=2*8, and the maximum stress-difference is 

'7 36X 2 # 8X2X1 5 = # 4 12X1 6 grammes per square centimeter. 

This stress-difference is, in British measure, 2*6 tons per square inch. 

If we suppose (as is not unreasonable) that it is 3 1 4 miles from crest to crest of the 
mountains, then the maximum stress will be reached at 50 miles below the surface. 

From Table VII. , § 9, it will be seen that if the materials of the earth at this depth 
of 50 miles had only as much tenacity as sheet lead, the mountain chains would sink 
down, but they would just be supported if the tenacity were equal to that of cast tin. 
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§ 7. On the stresses due to the even zonal harmonic inequalities. 

Having considered the two extreme cases where i is 2, and infinity, I pass now to 
the intermediate ones. As the odd zonal harmonics were not required for the investi- 
gation in the following section I have only worked out in detail the even ones. 

The surface of the sphere is now corrugated by a series of undulations parallel to 
the equator, and the altitude of the corrugations increases towards the poles. The 
form of the undulation in the neighbourhood of the equator is exhibited in Plate 19, 
fig. 3. 

The stress-difference is as before given by 

A=v / ^-R)M-4T 3 

To form this expression the series in (22) for R must be subtracted from the series 
for P. Since the (7 ? s and D's of Table II. have always the opposite signs from the A'& 
and jB's of Table I., this algebraic subtraction becomes a numerical addition of the 
numbers in these two tables. 

The results are given in the following table. 



Table IV. — The coefficients for expressing P-— R. 



* 


A —C 


A -Co 


+ 7-8431 


Aq~Cq 


^o-^o 


£«-£* 


Bt-Di 


■*e-»e 


2 


-2-2105 


-2-8421 




+ 2-5263 








4 

6 


-4-3137 


+ 2-7451 




+ 4-3922 


-7-5294 






-6-3636 


+ 42 


+ 35-6364 


-14-2545 


+ 6-4 


-48-8727 


+ 13-9636 


-21-5390 


8 


-8-3926 


+ 139-1411 


-38-8712 


-168-2061 


+ 8-4137 


-148-0806 


+ 188-4663 


10 
12 


-10-4115 


+ 318*2304 


-633-7449 


-478-0247 


+ 10-4252 


-329-2182 


+ 965-7051 


-491-6324 


-12-425 


+ 603-292 


-2623-009 


+ 552-212 


+ 12-434 


-616-315 


+ 3243-765 


-3806*018 



Then we have 

P-E=r*cotf l[(A -C ) + (A 2 -Cz) tan 3 1+ . . .] 

+aV-* cos- 2 l[(B -D Q ) + (Bs-D 2 ) tan 3 1+ . . . 

The materials for computing T have been already given in Table III. 

The series for P—-R and for 2T should now be squared and added together, but the 
result would be so complex that it is not worth while to proceed algebraically. 

At the equator T=0, and A=P — R, and P — R reduces to only two terms, 
whatever be the order of harmonic. 

By reference to (23) and (24) we see that at the equator 
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or 



it 



,i-2 



2(^ + 1)2 + 1 



'i(i + 2)(2i-l) 
i — 1 



a s — (i+l)(2i+3)r 



2 



ift + l)(2i + 3)aV-» 



(i + l)(2i + 3)-t 



r 



&, 



2 



+ ;: 



(i*-l)(2i + 3)_ 



(41) 



This value for A requires of course multiplication by appropriate factors involving 
the height of the continents and gravity. 

Even when i is no larger than 6. (41) 'differs but little from zV"" a (a 3 — r 2 ), at least for 
values of r not very nearly equal to a. 

A clearly reaches a maximum when 



W i-2 



\ a / 



i 



1+77 



O 

O 



(#»-l)(2i + 3)J 



For large values of i this maximum is nearly equal to 2{(i— 2)/i}**" 1 a* 
From these formulas the following results are easily obtained. 



Table V. (a). 



l-~ 


2 


4 


6 


8 


10 


12 


Maximiim value of A ... 
Value of rja when A is max. . 


2'526 



1-118 
714 


•959 
•819 


•894 
•867 


•859 
•895 


•836 
•913 



Plate 20, fig. 4, shows graphically the law of diminution of stress-difference for 
the even zonal harmonics, the vertical ordinates representing stress-difference and the 
horizontal ones the distances from the surface or from the centre of the globe. 

In order to find a numerical value of these maximum stress-differences which shall 
be intelligible according to ordinary mechanical ideas, I will suppose the height of 
each of the harmonics at the equator to be 1500 meters. On account of the small 
density of the superficial layers in the earth, this is very nearly the same as supposing 
that in the earth the maximum height of the continents above, and the maximum 
depth of the depressions below the mean level of the earth are each about 3000 meters. 
In the summary at the end I shall show that there is reason to believe that this is 
about the magnitude of terrestrial inequalities. 

Then by (33) we have to multiply the maximum stress -differences in the above 
table by 2(i-~l)ivh/(2i-\-l), in order to obtain the stress -differences for the supposed 
continents in grammes or tonnes per square centimeter. 

No"W*w==5 # 66 and A==l'5Xl0 5 according to the above hypothesis as to height of 
continent ; and the coefficient in i is of course different for each harmonic. 

By performing these multiplications I find the following results. 

MDCCOLXXXTT. 2 E 



210 



ME. G, H. DARWIN" ON THE STRESSES 



Table V. (6). — Maximum stress-differences due to harmonic continents and seas. 



Order of harmonic. 


2 


4 


6 


8 


10 


12 

•625 

3-96 
347 


Max. stress-difference, in metric tonnes per 
sq. cm. due to continents 1500 meters high 

Ditto in British tons per sq. inch, for same 
continents 


•858 

5-43 
f Centre 1 
\ of earth j 


•633 

4-01 
1146 


•626 

3-97 

725 


•625 

3-96 
532 


•625 

3-96 
420 


Depth in British miles at which this stress is 
attained 





N.B.—The continents referred to are supposed to he of the earth's mean density and are equivalent to actual 
continents of double the height. 

Thus far we have only considered the stress-differences at the equator immediately 
underneath the centres of the continents, but we must now see how they differ as the 
latitude of the place of observation increases. In order to attain this result a good 
deal of computation was necessary. 

For this purpose I calculated P— -R and 2T for a number of points and found the 
square root of the sum of these squares. As the computations were laborious, and as 
the results given in the following table are amply sufficient for the purpose in hand, 
I did not think it worth while to trace the changes to a greater depth than r=*7. 
Moreover the correctness of the last significant figures given cannot be guaranteed, 
although I believe that it is correct in most cases. 



Table VI. — Showing the stress-difference due to the several harmonic inequalities 

at various depths and in various latitudes. 



* 

•b 


Equator. 


Lat. 6°. 


Lat. 12°. 


• 

% 




Equator. 


Lat. 6°. 


Lat. 12°. 


fr=l. 


•316 


•316 


•316 




fr=l. 


♦021 


•015 


•000 


o J r='9 


•736 


•732 


•721 


8 


r=-9 


♦859 


•853 


•853 


1 


r='8 


1-112 


1108 


1-097 


r='8 


•798 


•795 


•797 




j— '7 


1-443 


1-440 


1-431 




jr='7 


•506 


•505 


•507 


fr=l. 
A, J r: =' 9 


•079 


♦074 


•061 




fr=l. 


•014 


•008 


•007 


•727 


•719 


♦700 


10 \ 


1 r=-9 


•857 


•854 


•860 


l 


r='8 


1-044 


1-038 


1-025 


r='S 


•631 


•630 


•635 




^=-7 


1-116 


1-113 


1-104 




^r='7 


•307 


•307 


•309 


fr=l. 
ft J r='9 


♦036 


♦031 


♦016 




fr=l. 


•010 


•003 


•019 


♦817 


•810 


•800 


12 < 


r=-9 


•827 


•824 


•835 


o s 


r=-8 


•953 


•949 


•945 


r='8 


•481 


•481 


•486 




w r=-7 


•788 


•786 


•785 




jr='7 


•179 


•179 


•181 



The numbers given in the column marked "equator" might be computed from (41). 
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and are those exhibited graphically in fig. 4 ; they are here given as a means of 
comparison with the numbers corresponding to latitudes 6° and 12°. 

The result to be deduced from this table is that the lines of equal stress-difference 
are very nearly parallel with the surface, and that it is for all practical purposes 
sufficient to know the stress-difference immediately under the centre of the continents. 

We have already seen in § 6 that for harmonics of infinitely high orders the lines of 
equal stress-difference are rigorously parallel with the mean surface. 



§ 8. On the stresses due to the iveight of an equatorial continent. 

The actual continents and seas on the earth's surface have not got quite the regular 
wavy character of the elevations and depressions which have been treated hitherto. 
The subject of the present section possesses therefore a peculiar interest for the 
purpose of application to the earth. Had I however foreseen, at the beginning, the 
direction which the results of this whole investigation would take, it is probable that 
I might not have carried out the long computations which were required for discussing 
the case of an isolated continent. But now that that end has been reached, it seems 
worth while to place the results on record. 

The function exp.[— cos 2 0/ sin 3 a] (where is colatitude) obviously represents an equa- 
torial belt of elevation, and I therefore chose it as the form of the required equatorial 
continent. This function has to be expanded in a series of zonal harmonics in order 
to apply the analytical solutions for the stresses produced by the weight of the 
continent. 

It is obvious by inspection that the odd zonal harmonics can take no part in the 
representation of the function, and it was on this account that I have above only 
worked out the cases of the even zonal harmonics. 

The multiplication of this function by the successive Legendre's functions, and 
integration over the surface of the sphere, are operations algebraically tedious, and 
wholly uninteresting, and I therefore simply give the results. 

I find then that if a =10°, and 

%= sin* 0-£ sin M 6 cos 3 + ^^ sin'-* cos 4 0-&c. 

Ci\ 41 

Then 



men 

2e-«- , »> ta --^ =i8 i5S3 +/8 4 s 4 +i8 fl s8+i8 8 S8+AoSio+A 8 5i2+ • • • 
where 

y8 =-3078, &=*3673, &=-3339, /3 6 =-2829,/3 8 =-2252, / 8 10 =-1688, 
/3 13 =-1193 



• ( 42 ) 



/3 is put on the left-hand side in order that the mean value of the function may he 
zero. The /3's obviously decrease very slowly, and as I stop with the 12th harmonic, 
the representation of the function is very imperfect. 

2 E 2 
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Plate 20, fig. 5, illustrates the results of the representation, the portion of a circle 
marked " mean level of earth " represents a meridional section of the earth ; the 
dotted curve marked " inequality to be represented " shows the true value of the function 
2exp.[— cos 2 6 1 sin 2 a]— /3 ; the curve marked " representation " shows the right-hand 
side of (42) stopping with the 12th harmonic; the second curve is the same without 
the 2nd harmonic constituent /3 3 %, and it is introduced for the reasons explained in the 
discussion and summary at the end. 

The equatorial value of the exponential function is 1792, that of the "representa- 
tion "is 1*497, and that of "the representation without the 2nd harmonic" is 1130. 

The polar value of the exponential is — '3078, that of the "representation" is 
— '0830, and that of "the representation without 2nd harmonic" is +"6516. This 
latter function thus gives us an equatorial continent and a pair of polar continents of 
less height. 

The extreme difference of height in the "representation" between the elevation at 
the equator and the depression at the pole is (r497 + *083)/& orl'58/i. I do not exactly 
know the extreme difference in the case where the 2nd harmonic is omitted, because 
I have not traced the inequality throughout, but it is probably about 1*3 or Vih 

Now let A,« be the numerical value, as computed for § 7, of the stress-difference due 
to the harmonic spheroid $, Then it has been shown above that the stress-difference 
due to the spheroid whose equation is r=a-\-h$i is —2(i—l)givhAi/(2i-{-l). 

Now stress-difference is a non-linear function of the component stresses P, R, T, and 
therefore the stress- difference due to a compound harmonic spheroid is not in general 
the sum of the stress-differences due to the constituent harmonic spheroids. At any 
point, however, where the principal stress-axes are all coincident in direction and where 
all the greater stress-axes coincide, and not a greater with a less, and where T=0, the 
stress- difference is linear and is the sum of the constituent stress-differences. This is 
the case at the equator for the present equatorial continent. 

Hence, if A be the stress-difference at the equator due to the spheroid, 

r=a+A(&s 2 +&s 4 + • . . +A2S12) 
We have 

A=-^[*&A s +fj8^^ . . . (43) 

In this formula the A/s are the numbers which were computed for drawing Plate 20, 
fig. 4, from the formula (41), namely 






r\* . 3 



1 — H +7^ 



a J ' (i 3 -l)(2?;+3)_ 



By using these computations I have drawn Plate 20, fig. 6. The vertical ordinates 
are — A+gwhy and the horizontal are the distances from surface or centre of the 
sphere, 
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The maxima in the two curves are merely found graphically, and the distances where 
the maxima are reached (viz. : 660 and 590 miles from the surface) are written down 
on the supposition that the radius of the sphere is 4000 miles. 

In the discussion in the second part of this paper, I have endeavoured to make an 
estimate of the proper elevation to attribute to these isolated continents ; so as to make 
the case, as nearly as may be, analogous to the earth. 

Although it appears impossible to make an accurate estimate, I conclude that it will 
not be excessive if we assume that the greatest difference of height, between the 
highest point in the equatorial elevation and the approximately spherical remainder of 
the globe, is 2000 meters. 

Accordingly for the representation I put 1*5 8A= 2000, and for the second curve 
l*4/i=2000 ; these give A=l*27 X 10 5 cm. and ^=1*4 X 10 5 cm. respectively. 

Taking w=5*66, then for the representation we have tvh— t 72 tonnes per square 
centimeter, and for the other curve ivh=*79 of the same units. The maximum stress- 
differences are '9lwh and *76wh respectively. 

Therefore for the equatorial table-land (called above the , representation) we have a 
maximum stress-difference of *66 metric tonnes per square cm. or 4*1 British tons per 
square inch ; and for the equatorial table-land balanced by a pair of polar continents 
(2nd harmonic omitted) we have a maximum stress -difference of '60 tonne per square 
cm. or 3*8 tons per square inch. 

I therefore conclude that our great continental plateaus produce a stress-difference of 
about 4 tons per square inch at a depth of 600 or 700 miles from the earth's surface. 



§ 9. On the strength of various substances. 

In order to have a proper comprehension of the strength which the earth's mass 
must possess in order to resist the tendency to rupture, produced by the unequal 
distribution of weights on the surface, it is^ necessary to consider the results of 
experiments. 

Rankine* gives a large number of results obtained by various experimenters, and 
Sir William Thomson also gives similar tables in his article on ' Elasticity/ 1 

Amongst other constants Sir William gives Young's modulus and the greatest 
elastic extension. If the materials of a wire remain perfectly elastic when the wire is 
just on the point of breaking under tension, then the product of Young's modulus 
into the greatest elastic extension should be equal to what is called the tenacity, 
which is defined as the breaking tension per square centimeter of the area of the wire. 

* ' Useful Rules and Tables:' Griffin, London, 1873, p. 191, et seq. 

t ' Elasticity : ' Black, Edinburgh, 1878. This is the article from the ' Encyclopedia Britannica,' but 
it is also published as a separate wwk. 
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If however a permanent set begins to take place before the wire breaks, this product 
should be less than the tenacity. I do not see how it can ever be greater, unless 
there be a marked departure from Hook's law "ut tensio sic vis;" or different sets of 
experiments with the same class of material might make it seem greater. In some of 
the results given by Sir William Thomson the product of modulus and elastic 
extension is however greater than tenacity. 

Ordinary experience would lead one to suppose that such materials as lead and 
copper would undergo a considerable stress beyond the limits of perfect elasticity, 
before breaking. It is surprising therefore to see how nearly identical this product is 
to the tenacity — indeed in the case of lead absolutely identical, as may be seen in the 
table below. 

With regard to the earth we require to know what is the limiting stress-difference 
under which a material takes permanent set or begins to flow, rather than the stress- 
difference under which it breaks ; for if the materials of the earth were to begin to 
flow, the continents would sink down and the sea bottoms rise up. 

It will be seen from the definition of tenacity given above that it is the rupturing 
stress- difference for tensional stresses. There is no word specially applied to rupturing 
stress- difference under pressure. 

I am inclined to think that for the purposes of this investigation these tables in 
most cases rate the strength of materials somewhat too highly; for it seems probable 
that a permanent set would be taken, if a material were subjected for a long time to a 
stress-difference, which is a considerable fraction of the limiting value. We are likely 
to know more on this point in some years time when the wires hung by Sir William 
Thomson in the tower of Glasgow University have been subjected to several years of 
tension. However this may be I give the results of some of the experiments as 
collected and quoted by Sir William Thomson and the late Professor Kankine. 
The first table of tenacity, except the results denoted by the letter R, are taken from 
Sir William Thomson. The second table of crushing stress- difference is taken 
entirely from Rankine. The multiplications and reductions to different units I have 
done myself. 
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Table VII. — Limiting stress-difference. 



Produced by tension. 



Material. 



Sheet lead . 
Cast tin 



J5 



Wood (ash) 
Cast brass . 



5> 



iron 



Drawn copper 
English steel piano 

forte wire . . 
(R) Brick, cement 
(R) Glass . . . 
(R) Slate . . . 



Breaking 
stress- 
difference in 
metric tonnes 
per square 
centimeter. 



•23 

•416 

•325 (R) 
1-20 
1-27 
•94 to 2-04 
4-10 

23-62 
020 to -021 

'66 
•68 to -90 



Stress-difference at which 
permanent set begins in — 



Metric tonnes 
per square 
centimeter. 



British tons 

per square 

inch. 



•23 
•417 


1-20 
1-27 

•14 to 1-87 
4-00 


23-56 

• • 


• • 

• • 



1-46 
2-64 
2-06 (R) 
7-61 
8-05 
7-23 to 11- 
25-36 

149-6 
•125 to -134 

4-20 
4-3 to 5-7 



Produced by crushing. 



Material. 



Strong red brick 
Strong sandstone 
(F) Strong lime 

stone . . . 
Marble 
Granite . . 
(F) Granite (Mount 

Sorrel) . . 
(F) Grauwacke . 
Ash (along the grain) 
Cast brass . . 
Wrought iron . 



Breaking stress-difference 
in — 



Metric tonnes 
per square 
centimeter. 



•077 
•39 

•60 

•39 

•39 to 



77 



•905 
1-19 
•63 
•73 

2-52 to 2-84 



British tons 

per square 

inch. 



•49 
2-45 

3-80 

2-45 

2-45 to 4-91 

5-74 
7-54 
4-02 
4-60 
16 to 18 



Note. — The second and third columns give the product of Young's modulus into greatest elastic exten- 
sion, and this should give the stress-difference when permanent set begins. Rankine does not give the 
data for this quantity, but the breaking stress- difference is given in both metric and British units, the 
latter being in the third column. 

In the second half of the table the results marked F are from Sir William Fairbairn's experiments. 

The only cases in these two tables in which we have the opportunity of comparing 
the strength for resisting the stress-difference, when produced in the two manners, is 
for the materials cast brass and ash; in both cases we see that the substance is 
considerably weaker for crushing than for tension. 

I should be inclined to suppose that the crushing strength is more nearly the datum 
we require for the case of the stresses in the earth. 

In the first half of the table we probably see the effect of permanent set in the 
cases of copper and pianoforte wire (compare 4*00 with 4*10, and 23*56 with 23*62), 
but it is surprising that the contrast between the two columns is not more marked. 



§ 10. On the case when the elastic solid is compressible. 

It appears desirable to know how far the results of the preceding investigation may 
differ, if the elastic solid be compressible. According to the views of Dr. Hitter, 
referred to in the summary, this must be largely the case. 

As I did not examine this point until after the above was completed, it seemed 
preferable to make a fresh beginning, rather than to modify the whole investigation. 
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The process is however so closely analogous, that it presents no difficulty, and may be 
dismissed shortly. I shall accordingly follow closely the processes of §§ 1 and 2. 

If a solid be very compressible it takes a comparatively small hydrostatic pressure 
to produce a given amount of compression : that is to say, although the compressibility 
is large, the modulus of compressibility is small compared with that of rigidity. The 
modulus of compressibility I shall call the incompressibility. In the preceding investi- 
gation the converse was the case, for the incompressibility was infinite compared with 
the rigidity. By the definitions of co and v in § 1, the incompressibility 



k= 



O) 



-v. 



It will be found convenient to use h and co as the two moduli, which define the 
nature of the elastic solid. 

In the denominators in E h F h Gi of (1), the expression (2(i+l) 3 + !)<*>"-( 2i+l)v 

occurs, this I shall call K, in analogy with i". 

Then 

K=2i(i-~l)G)+ 3(2i+l)k 

If we develop the last differential coefficient in the expression (1) for a, we find 



2vKa=3k {- 



% 



i-1 



■oft — T * 



0X dW , . 



dx 



(a i -r A )^-{-2xW 



dx 



Also 



KB=iW . 



. * « » « 



• • • ( 44 ) 



* * • . » 1 ^t%J J 



and 



P=(3£-2co)S+2i; 

t_ ( clrx i d ^ 

\dz dx j 



da 



d 



x 



r 



6 • » 



(46) 



Differentiating (44) with regard to x, and substituting in (46) we find 



JLx. X. — — O fC 



% 



\d?W rt dW , , TT " 
--. a *_ r - )— — — 2w~ T --\-iW 

.i— 1 / oar ao? 



r2. 






(47) 



Again, differentiating (44) with regard to % 9 writing down dy/dx by symmetry, and 
adding the two together, we have from (46) 



Jtl JL — d/C 



% 



i — \ 



9 AcPW dW dW 
or— r-'hr-; %-r- — x~ 



Jdxdz 



dx 



dz 



+^(a 3 —-r 2 ) 



Cl/iX/CvZ 



(48) 



From (47) a,nd (48) the other stresses may be written down by cyclic changes of 
letters. 

Now let us suppose that the incompressibility is small compared with the rigidity. 
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Then o) becomes nearly equal to ^v, and h is small compared with v and co. Also K 
becomes nearly equal to 2i(i—i)<a. Hence in this case 



1 cP W 

2(i-iy l } da? 

1 cPW 

2(i-l^ a ~ r 'dxdz 



(49) 



Now if we suppose W to be a zonal harmonic of degree i, and only consider the 
state of stress at the equator, immediately underneath the centre of the elevation, 
then T is zero and by (15) 



R=-- 7 Arr(a 3 -r s y'" 3 

2(1 — l) x ' 



iiDCl 



A=P-R=^-^(a s -^y- 3 (50) 

If i be large A=i(a s — r^)r*" 2 . 

In §7 we have found this identical result, under the like conditions, when the solid 
is incompressible. 

Now take the case of the 2nd harmonic, so that 



And we have 



W=r\\- cos 2 0)=i(a*+y*-2z*) 

P= -|(a 2 — r 2 ) 
Q= i(a 2 — r 2 ) 
K=-f(a 2 -r 2 ) 
T = 

Thus throughout the spheroid, the principal stress-axes are parallel and perpen- 
dicular to the polar axis ; also 

A = (a 3 -r 3 ) 

Hence the central stress-difference is a 2 . If the solid be incompressible it was found 
to be 8a 2 . Hence infinite compressibility largely relieves the stress-difference due to 
ellipticity of figure. Next take the case when the harmonic is of the 4th order. 

Then at the equator we have by (50) 



6 



r% (a? — r 2 ) = ^r 2 (a 2 — r 2 ) 



The maximum is reached, when r % =\a? or r=*707a, and is equal to |= 11*67. 
Comparing this with the Table V. (a) § 7, we see that infinite compressibility makes 
very little difference, for 707 differs little from 714 and 1167 from 1118. 
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It will now be shown that whatever may be the compressibility of the solid, we get 
identically the same solution in the case when the harmonic is of an infinitely high 
order. This is the problem of the harmonic mountains and valleys corrugating a mean 
level surface, which was considered in § 6, The same notation will be adopted here. 

Both i and a are infinite, x becomes a— £ and ia 2 /(i—l)—r & =a?—r*=2a£;. If the 
substitutions here suggested be made in (47) and (48), it will be found that the terms 
with (o as a factor are multiplied by ai (two infinites), whilst none of the terms with k 
as a factor involve more than one infinite. Hence the latter terms are negligeable 
compared with the former. 

Again i being infinite, K=2i 2 a). Thus if i and a be infinite (47) and (48) reduce to 

1 d?W 1 d*W 1 d?W 

P= 2T^ 2a ^> R= ^' 2a ^' T =-2^" 2a ^fe 



but a/i=b, therefore 



P =«5F- *=*?-£• T =-<7fi 



In (36) we found the effective potential, which produced the same state of stress as 
the weight of the mountains whose equation was -~g=zh cos z/b ; the result was 

W= —gioh e~~^' b cos ~ 
Hence in the present case 



* 'V 



P— : — gwhj-e~V b cos'j, R= givhye"^ cos ~, Ti=giuhje'~ l/b sin- 

These results are identical with equations (38), which gave the stresses when the 
elastic solid was incompressible. 

It may be concluded from this that, except for the case of the 2nd harmonic 
inequality, compressibility makes very little difference, and for the higher harmonics it 
makes no difference at all. 

II. 

SUMMARY AND DISCUSSION. 

The existence of dry land proves that the earth's surface is not a figure of 
equilibrium appropriate for the diurnal rotation. 

Hence the interior of the earth must be in a state of stress, and as the land does 
not sink in, nor the sea-bed rise up, the materials of which the earth is made must be 
strong enough to bear this stress. 

We are thus led to inquire how the stresses are distributed in the earth's mass, and 
what are their magnitudes. These points cannot be discussed without an hypothesis 
as to the interior constitution of the earth. 
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In this paper I have solved a problem of the kind indicated for the case of a 
homogeneous incompressible elastic sphere, and have applied the results to the case 
of the earth. 

It may of course be urged that the earth is not such as this treatment postulates. 

The view which was formerly generally held was that the earth consists of a solid 
crust floating on a molten nucleus. It has also been lately maintained by Dr. August 
Hitter in a series of interesting papers that the interior of the earth is gaseous/* 
A third opinion, contended for by Sir William Thomson, and of which I am myself 
an adherent, is that the earth is throughout a solid of great rigidity ; he explains 
the flow of lava from volcanoes either by the existence of liquid vesicles in the interior, 
or by the melting of solid matter, existing at high temperature and pressure, at points 
where diminution of pressure occurs. 

There is another consideration, which is consistent with Sir William Thomson's 
view, and which was pointed out to me by Professor Stokes. It may be that undei;- 
neath each continent there is a region of deficient density; then underneath this 
region there would be no excess of pressure. 

For the present investigation it is to some extent a, matter of indifference as to 
which of these views is correct, for if it is only the crust of the earth which possesses 
rigidity, or if Professor Stokes's suggestion of the regions of deficient density be 
correct, then the stresses in the crust or in the parts near the surface must be greater 
than those here computed — enormously greater if the crust be thin,t or if the region 
of deficient density be of no great thickness. 

With regard to the property of incompressibility which is here attributed to the 
elastic sphere, it appears from §10 that even if we suppose the elastic solid to be very 
highly compressible, yet the results with regard to the internal stresses are almost the 
same as though it were incompressible. I think tlie hypothesis of great incompressi- 
bility is likely to be much nearer to the truth than is that of great compressibility. 
I shall therefore adhere to the supposition of infinite incompressibility, bearing in 
mind that even great compressibility would not much affect most of the results. 

I take then a homogeneous incompressible elastic sphere, and suppose it to have the 

* 'Anwendung der mechanischen Warmetheorie auf kosmologische Probleme.' Carl Rumpler, 
Hannover, 1879. This is a reprint of six papers in Wiedemann's Annalen. 

Dr. Ritter contends that the temperature in the interior of the planet is above the critical tempera- 
ture and that of dissociation for all the constituents, so that they can only exist as gas. Data are wanting 
with regard to the mechanical properties of matter at, say 10,000° Fahr., and a pressure of many tons to 
the square inch. Is it not possible that such "gas" may have the density of mercury and the rigidity 
and tenacity of granite ? Although such a conjectural " gaseous " solid might possess high rigidity, it 
almost certainly would have great compressibility : but it is proved in § 10 that the compressibility will 
make exceedingly little difference in the result of the present investigation excepting in the case of the 
2nd harmonic inequality. 

t The evaluation of the stresses in a crust, with fluid beneath, would be tedious, but not more difficult 
than the present investigation. I may perhaps undertake this at some future time. 

2 F 2 
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power of gravitation and to be superficially corrugated. In consequence of mathe- 
matical difficulties the problem is here only solved for the particular class of surface 
inequalities called zonal harmonics, the nature of which will be explained below. 

Before discussing the state of stress produced by these inequalities, it will be con- 
venient to explain the proper mode of estimating the strength of an elastic solid under 
stress. 

At any point in the interior of a stressed elastic solid there are three lines mutually 
at right angles, which are called the principal stress-axes. Inside the solid at the 
point in question imagine a small plane (say a square centimeter or inch) drawn 
perpendicular to one of the stress-axes ; such a small plane will be called an inter- 
face/* The matter on one side of the ideal inter-face might be removed without 
disturbing the equilibrium of the elastic solid, provided some proper force be applied 
to the inter-face ; in other words, the matter on one side of an inter-face exerts a force 
on the- matter on the other side. Now a stress-axis has the property that this force 
is parallel to the stress-axis to which the inter-face is perpendicular. Thus along a 
stress-axis the internal force is either purely a traction or purely a pressure. Treating 
pressures as negative tractions, we may say that at any point of a stressed elastic 
solid, there are three mutually perpendicular directions along which the stresses are 
purely tractional. The traction which must be applied to an inter-face of a square 
centimeter in area, in order to maintain equilibrium when the matter on one side of 
the inter-face is removed, is called a principal stress, and is of course to be measured 
by grammes weight per square centimeter. 

If the three stresses be equal and negative, the matter at the point in question is 
simply squeezed by hydrostatic pressure, and it is not likely that in a homogeneous 
solid any simple hydrostatic pressure, absolutely equal in all directions, would ever 
rupture the solid. The effect of the equality of the three stresses when they are 
positive and tractional is obscure, but at least physicists do not in general suppose 
that this is the cause of rupture when a solid breaks. 

If the three principal stresses be unequal, one must of course be greatest and one 
least, and there is reason to suppose that tendency of the solid to rupture is to be 
measured by the difference between these principal stresses. 

In one very simple case we know that this is so, for if we imagine a square bar, of 
which the section is a square centimeter, to be submitted to simple longitudinal ten- 
sion, then two of the principal stresses are zero (namely, the stresses perpendicular to 
the faces of the rod), and the third is equal to the longitudinal traction. The traction 
under which the rod breaks is a measure of its strength, and this is equal to the 
difference of principal stresses. 

If at the same time the rod were subjected to great hydrostatic pressure, the break- 
ing load would be very little, if at all affected ; now the hydrostatic pressure subtracts 

* This term is due to Professor James Thomson. 
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the same quantity from all three principal stresses, but leaves the difference between 
the greatest and least principal stresses the same as before. 

Difference of principal stresses may also be produced by crushing. 

In this paper I call the difference between the greatest and least principal stresses 
the " stress-differ ence," and I say that, if calculation shows that the weight of a certain 
inequality on the surface of the earth will produce such and such stress-difference at such 
and such a place, then the matter at that place must be at least as strong as matter 
which will break when an equal stress-difference is produced by traction or crushing. 

I shall usually estimate stress-difference by metric tonnes (a million grammes) per 
square centimeter, or by British tons per square inch. 

In Table VII., § 9, are given the experimentally determined values of the breaking 
stress-difference for various substances. The table is divided into two parts, in the 
former of which the stress-difference was produced by tension, and in the latter by 
crushing. It is not necessary here to advert to the difference in meaning of the 
numbers given in the first column and those given in the two latter columns in the 
first half of the table. 

The cases of wood and cast brass are the only ones where a comparison is possible 
between the two. breaking stress-differences, as differently produced. It will be seen 
that the material is weaker for crushing than for tension. For the reasons given in 
that section, I am inclined to think that these tables rate the strength of the materials 
somewhat too highly for the purposes of this investigation. I conceive that the 
results derived from crushing are more appropriate for the present purpose than those 
derived from tension ; and fortunately the results for various kinds of rocks seem to 
have been principally derived from crushing stresses. 

This table will serve as a means of comparison with the numerical results derived 
below, so that we shall see, for example, whether or not at 500 miles from the surface 
the materials of the earth are as strong as granite. 

We may now pass to the mathematical investigation. It appears therefrom that 
the distribution of stress-difference is quite independent of the absolute heights and 
depths of the inequalities. Although the questions of distribution and magnitude of 
the stresses are thus independent, it will in general be convenient to discuss them 
more or less simultaneously. 

The problem has only been solved for the class of superficial inequalities called 
zonal harmonics, and their nature will now be explained. 

A zonal harmonic consists of a series of undulations corrugating the surface in 
parallels of latitude with reference to some equator on the globe ; the number of the 
undulations is estimated by the order of the harmonic. The harmonic of the second 
order is the most fundamental kind, and consists of a single undulation forming an 
elevation round the equator, and a pair of depressions at the poles of that equator ; it 
may also be defined as an elliptic spheroid of revolution, and the absolute magnitude 
is measured by the ellipticity of the spheroid. 
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If the order of the harmonic be high, say 30 or 40, we have a regular series of 
mountain chains and intervening valleys running round the sphere in parallels of 
latitude. 

For the sake of convenience I shall alwaj^s speak as though the equator were a 
region of elevation, but the only effect of changing elevations into depressions, and 
vice versd, is to diametrically reverse the directions of all the stresses. 

The harmonics of the orders 2, 6, 10, &a, have depressions at the poles of the 
sphere; those of orders 4, 8, 12, &c>, have elevations at the pole. 

The harmonic of the fourth order consists of an equatorial continent and a pair of 
circular polar continents, with an intervening depression. That of the sixth order 
consists of an equatorial continent and a pair of annular continents in latitudes 
(about) 60° on one and the other side of the equator. The 8th harmonic brings 
down these new annular continents to about latitude 45°, and adds a pair of polar 
continents ; and so on. 

By a continuation of this process the transition to the mountain chains and valleys 
is obvious. 

In § 5 the case of the 2nd harmonic is considered. As above explained the 
sphere is deformed into a spheroid of revolution. The investigation also applies to 
the case of a rotating spheroid, such as the earth, with either more or less oblateness 
than is appropriate for the figure of equilibrium. 

The lines throughout a meridional section of the spheroid along which the stress- 
difference is constant are shown in Plate 19, fig. 1, and the numbers written on the 
curves give the relative magnitude of the stress- difference. 

It is remarkable that the stress-difference is the same all over the surface. In the 
polar regions the stress-difference diminishes as we descend into the spheroid and then 
increases again ; in the equatorial regions it always increases as we descend. The 
maximum value is at the centre, and there the stress-difference is eight times as great 
as at the surface. 

If the elastic solid be highly compressible the stress-differences are not nearly so 
great as on the hypothesis of incompressibility. In all the other cases considered in 
this paper compressibility makes practically no difference in the results. 

On evaluating the stress-difference, on the hypothesis of incompressibility, arising 
from given ellipticity in a spheroid of the size and density of the earth, it appears 
that if the excess or defect of ellipticity above or below the equilibrium value 
(namely g-g-g for the homogeneous earth) were 7^0 q-, then the stress-difference at the 
centre would be 8 tons per square inch, and accordingly, if the sphere were made of, 
material as strong as brass (see Table VII.), it would be just on the point of rupture. 
Again if the homogeneous earth, with ellipticity g-g-g, were to stop rotating, the 
central stress- difference would be 33 tons per square inch, and it would rupture if 
made of any material excepting the finest steel. 
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A rough calculation 4 ' will show that if the planet Mars has ellipticity -^- - (about 
twice the ellipticity on the hypothesis of homogeneity) the central stress-difference 
must be 6 tons per square inch. It was formerly supposed that the ellipticity of the 
planet was even greater than g^j, an( i even if the latest telescopic evidence had not 
been adverse to such a conclusion, we should feel bound to regard such supposed 
ellipticity with the greatest suspicion, in the face of the result just stated. 

The state of internal stress of an elastic sphere under tide-generating forces is 
identical with that caused by ellipticity of figure. t Hence the investigation of § 5 
gives the distribution of stress-difference caused in the earth by the moon's attraction. 
In Plate 1 9, fig. 1, the point called " the pole" is the point where the moon is in the zenith. 

Computation shows that the stress-difference at the surface, due to the lunar ti de- 
generating forces, is 16 grammes per square centimeter, and at the centre eight times 
as much. These stresses are considerable, although very small compared with those 
due to terrestrial inequalities, as will appear below. 

In § 6 the stresses produced by harmonic inequalities of high orders are considered. 
This is in effect the case of a series of parallel mountains and valleys, corrugating a 
mean level surface with an infinite series of parallel ridges and furrows. In this case 
compressibility makes absolutely no difference in the result, as shown in § 1 0. 

It is found that the stress-difference depends only on the depth below the mean 
surface, and is independent of the position of the point considered with regard to 
ridge and furrow ; the direction of the stresses does however depend on this latter 
consideration. 

In Plate 19, fig. 2, is shown the law by which the stress-difference increases and then 
diminishes as we go below the surface. The vertical ordinates of the curve indicate the 
relative magnitude of the stress-difference, and the horizontal ones the depth below the 
surface. The depth OL on the figure is equal to the distance between contiguous 
ridges, and the figure shows that the stress-difference is greatest at a depth equal to 

The greatest stress-difference depends merely on the height and density of the. 
mountains, and the depth at which it is reached merely on the distance from ridge to 

ridge. 

Numerical calculation shows that if we suppose a series of mountains, whose crests 
are 4000 meters or about 13,000 feet above the intermediate valley-bottoms, formed 
of rock of specific gravity 2*8, then the maximum stress-difference is 2*6 tons per 
square inch (about the tenacity of cast tin) ; also if the mountain chains are 314 miles 
apart the maximum stress-difference is reached at 50 miles below the mean surface. 

* The data for the calculation are : Ratio of terrestrial radius to Martian radius 1*878. Ratio of 
Martian mass to terrestrial mass '1020. Whence ratio of Martian gravity to terrestrial gravity '3596. 
Central stress-difference, due to ellipticity e, 996e tons per square inch. " Homogeneous " ellipticity of 
Mars T |~e ; and fff equal to 6. 

f This is subject to certain qualincations noticed in § 5. 
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It may be necessary to warn the geologist that this investigation is approximate in a 
certain sense, for the results do not give the state of stress actually within the mountain 
prominences or near the surface in the valley -bottoms. The solution will however be 
very nearly accurate at some five or six miles below the valley-bottoms. The solution 
shows that the stress-difference is nil at the mean surface, but it is obvious that both 
the mountain masses and the valley-bottoms are in some state of stress. 

The mathematician will easily see that this imperfection arises, because the problem 
really treated is that of an infinite elastic plane, subjected to simple harmonic tractions 
and pressures. 

To find the state of stress actually within the mountain masses would probably be 
difficult. 

The maximum stress-difference just found for the mountains and valleys obviously 
cannot be so great as that at the base of a vertical column of this rock, which has a 
section of a square inch and is 4000 meters high. The weight of such a column is 7*1 
tons, and therefore the stress- difference at the base would be 7*1 tons per square inch. 
The maximum stress-difference computed above is 2*6, which is about three-eighths 
of 7*1 tons per square inch. Thus the support of the contiguous masses of rock, in the 
case just considered, serves as a relief to the rock to the extent of about five-eighths 
of the greatest possible stress-difference. This computation also gives a rough estimate 
of the stress-differences which must exist if the crust of the earth be thin. It is 
shown below that there is reason to suppose that the height from the crest to the bottom 
of the depression in such large undulations as those formed by Africa and America is about 
6000 meters. The weight of a similar column 6000 meters high is nearly 11 tons. 

In §71 take the cases of the even zonal harmonics from the 2nd to the 12th, but 
for all except the 2nd harmonic only the equatorial region of the sphere is considered. 

Plate 19, fig. 3, shows an exaggerated outline of the equatorial portion of the inequali- 
ties ; it only extends far enough to show half of the most southerly depression, even for 
the 12th harmonic. It did not seem worth while to trace the surfaces of equal stress- 
difference throughout the spheroid, but the laborious computations are carried far 
enough to show that these surfaces must be approximately parallel to the surface of 
the mean sphere. It is accordingly sufficient to find the law for the variation of stress-. 
difference immediately underneath the equatorial belt of elevation. It requires com- 
paratively little computation to obtain the results numerically, and the results of the 
computation are exhibited graphically in Plate 20, fig. 4. 

Table V. (b), § 7, gives the maximum stress-differences, resulting from these several 
inequalities, computed under conditions adequately noted in the table itself. It will 
be convenient to postpone the discussion of the results. 

In § 8 I build up out of these six harmonics an isolated equatorial continent. The 
nature of the elevation is exhibited in Plate 20, fig. 5, in the curve marked " represen- 
tation;" no notice need be now taken of the dotted curve. This curve exhibits a belt 
of elevation of about 1 5° of latitude in semi-breadth, and the rest of the spheroid is 
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approximately spherical. This kind of elevation requires the 2nd as one of its harmonic 
constituents, and this harmonic means ellipticity of the whole globe. Now it may 
perhaps be fairly contended that on the earth we have no such continent as would 
require a perceptible 2nd harmonic constituent. I therefore give in Plate 20, fig. 5, a 
second curve w^hich represents an equatorial belt of elevation counterbalanced by a pair 
of polar continents in such a manner that there is no second harmonic constituent. 

I have not attempted to trace the curves of equal stress-difference arising from 
these two kinds of elevation, but I believe that they will consist of a series of much 
elongated ovals, whose longer sides are approximately parallel with the surface of the 
globe, drawn about the maximum point in the interior of the sphere at the equator. 
The surfaces of equal stress-difference in the solid figure will thus be a number of 
flattened tubular surfaces one within the other. 

"At the equator however the law of variation of stress-difference is easy to evaluate, 
and Plate 20, fig. 6, shows the results graphically, the vertical ordinates representing 
stress-difference and the horizontal the depths below the surface. The upper curve in 
Plate 20, fig. 6, corresponds with the "representation curve" of Plate 20, fig. 5, and 
the lower curve with the case where there is no 2nd harmonic constituent. 

The central stress-difference, which may be observed in the upper curve, results 
entirely from the presence of the 2nd harmonic constituent in the corresponding 
equatorial belt of elevation. 

The maximum stress-differences in these two cases occur at about 660 and 590 miles 
from the surface respectively. 

We now come to perhaps the most difficult question with regard to the whole 
subject — namely, how to apply these results most justly to the case of the earth. 

The question to a great extent turns on the magnitude and extent of the superficial 
inequalities in the earth. As the investigation deals with the larger inequalities, it 
will be proper to suppose the more accentuated features of ridges, peaks, and holes to 
be smoothed out. 

The stresses caused in the earth by deficiency of matter over the sea beds are the 
same as though the seas were replaced by a layer of rock, having everywhere a 
thickness of about -|^|- or nearly -^ of the actual depths of sea. 

The surface being partially smoothed and dried in this manner, we require to find 
an ellipsoid of revolution which shall intersect the corrugations in such a manner that 
the total volume above it shall be equal to the total volume below it. 

Such a spheroid may be assumed to be the figure of equilibrium appropriate to the 
earth's diurnal rotation ; if it departs from the equilibrium form by even a little, then 
we shall much underestimate the stress in the earth's interior by supposing it to be a 
form of equilibrium. 

Professor Bruns has introduced the term " geoid " to express any one of the " level " 
surfaces in the neighbourhood of the earth's surface, and he endeavours to form an 
estimate of the departure of the continental masses and sea-bottoms from some mean 
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geoid/* From the geodesic point of view the conception is valuable, but such an 
estimate is scarcely what we require in the present case. The mean geoid itself will 
necessarily partake of the contortions of the solid earth's surface, even apart from 
disturbances caused by local inequalities of density, and thus it cannot be a figure of 
equilibrium. 

Thus, even if we were to suppose that the solid earth were everywhere coincident 
with a geoid — which is far from being the case — a state of stress would still be 
produced in the interior of the earth. 

An example of this sort of consideration is afforded by the geodesic results arrived 
at by Colonel Clarke, R.E.,t who finds that the ellipsoid which best satisfies 
geodesic measurement, has three unequal axes, and that one equatorial semi-axis is 
1524 feet longer than the other. Now such an ellipsoid as this, although not exactly 
one of Bruns' geoids, must be more nearly so than any spheroid of revolution ; and 
yet this inequality (if really existent, and Colonel Clarke's own words do not express 
any very great confidence) must produce stress in the earth. Colonel Clarke's 
results show an ellipticity of the equator equal to t"3T3T> an d this in the homogeneous 
elastic earth will be about equivalent to ellipticity fy^oo J suc h ellipticity would 
produce a central stress-difference of -^Tooo or nearly one-third of a British ton per 
square inch. 

From this discussion it may, I think, be fairly concluded that if we assume the sea- 
level as being the figure of equilibrium and estimate the departures therefrom, we shall 
be well within the mark. 

The average height of the continents is about 350 meters (1150 feet), and the 
average depth of the great oceans is in round numbers 5000 meters (16,000 feet) ; 
but the latter datum is open to much uncertainty.! When the sea is solidified into 
rock the 5000 meters of depth is reduced to 3200 meters below the actual sea-level. 
Thus the average effective depression of sea-bed is about nine times as great as the 
average height of the land. I shall take it as exactly nine times as great, and put the 
depth as 3150 meters; but it is of course to be admitted that perhaps eight and 
perhaps ten might be more correct factors. 

In the analytical investigation . of this paper the outlines of the vertical section of the 
continents and depressions are always sweeping curves of the harmonic type, and the 
magnitude of the elevations and depressions are estimated by the greatest heights 
and depths, measured from a mean surface which equally divides the two. 

We have already supposed the outlines of continents and sea-beds to have been 
smoothed down into sweeping curves, which we may take as being, roughly speaking, 
of the harmonic type. The smoothing will have left the averages unaffected. 

* 'Die Eigur der Erde.' Yon Dr. H. Bruns, Berlin: Stankiewicz, 1878. 
t Phil. Mag., Aug., 1878. 

J In a previous paper, " Geological Changes, &c.," Phil. Trans., Yol. 167, Part L, p. 295, I have 
endeavoured to discuss this subject, and references to a few authorities will be found there. 
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The averages are not however estimated from a mean spheroidal surface, but from 
one which is far distant from the mean. 

The questions now to be determined are as follows : — What is the proper greatest 
height and depression, estimated from a mean spheroid, which will bring out the above 
averages estimated from present sea-level, and what is the position of the mean 
spheroid with reference to the S ea-level. ^ 

From the solution of the problem considered in the note below, 41 it appears that, if 

* Conceive a series of straight harmonic undulations corrugating a mean horizontal surface, and 
suppose them to be flooded with water. This will represent fairly well the undulations on the dried earth, 
and the water-level will represent the sea-level. 

Suppose that the average heights and depths of the parts above and below water are known, and 
that it is required to find the position of the mean horizontal surface with reference to the water-level, 
and the height of the undulations measured from that mean surface. 

Take an origin of coordinates in the water-level, the axis of x in the water-level and perpendicular to 
the undulations, and the axis of y measured upwards. 

Let 

y=h(cos x~ cos a) 
be the equation to the undulations. 



1 



* +a , h 



The average height of the dry parts is clearly 3- ydx or -(sin a— a cos a). Similarly the average 

h h 

depth below water is ----[sin (tt— a)~ (ir—a) cos (tt— cl)~] or -^- sin a+ (ir—a) cos a 

If the latter average be p times as great as the former 

ph cos a(- tan a— 1 )=h cos a( tan a + I ) 

\a / V— a ) 

This is an equation for determining a. 

Now I find that a=34° 30' gives p= 8*983, which corresponds very nearly with p=9 of the text above. 
This value of a corresponds with an average equal to '1165& for the height above water, and 1'0469/i 
for the depth below water. Now if we put 

1-0469^=3150 meters 
which gives *11657i=350 meters very nearly, 
we have ^=3009 meters. 

The depth below water-level of the mean level is h cos 34° 30' or 2480 meters. 

The greatest height of the dry part above the water-level is 3009—2480 or 429 meters, and the 
greatest depth of the submerged part below water-level is 3009 + 2480 or 5489 meters. 

[After the proof-sheets of this paper had been corrected, Professor Stokes pointed out to me that, 
according to Rigaud (Cam. Phil. Soc, vol. 6), the area of land is about four-fifteenths of the whole area 
of the earth's surface. Now, in the ideal undulations we are here considering the area above water is 
about one-tenth of the whole area ; hence in this respect the analogy is not satisfactory between these 
undulations and the terrestrial continents. If I have not considerably over-estimated the average depth 
of the sea (and I do not think that I have done so), the discrepancy must arise from the fact that actual 
continents and sea-beds do not present in section curves which conform to the harmonic type ; there must 
also be a difference between corrugated spherical and plane surfaces. 

The geological denudation of the land must, to some extent, render our continents flat-topped. — Added 
May Mh, 1882.] 
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the continents and sea-beds have sections which are harmonic curves, then if we 
take, — 

The mean level bisecting elevations and depressions as 2480 meters (8150 feet) 
below the sea-level, and the greatest elevation and depression from that mean level as 
3009 meters (9840 feet), it results that the average height of the land above sea-level 
is 350 meters and the average depression of dried sea-bed is 3150 meters. 

It thus appears that 3000 meters would be a proper greatest elevation and de- 
pression to assume for the harmonic analysis of this paper, if the earth were 
homogeneous. But as the density of superficial rocks is only a half of the mean 
density of the earth, I shall take 1500 meters as the greatest elevation and depression 
from the mean equilibrium spheroid of revolution. 

- It is proper here to note that the height of the undulations of elevation and depres- 
sion in the zonal harmonic inequalities is considerably greater towards the poles than 
it is about the equator; it might therefore be maintained that by making 1500 meters 
the equatorial height, we are taking too high an estimate. But the state of stress 
caused in the sphere at any point depends very much more on the height of the 
inequality in the neighbourhood of a superficial point immediately over the point 
considered, than it does on the inequalities in remote parts of the sphere. 

Now in all the inequalities, except the 2nd harmonic, I have considered the state of 
stress in the equatorial region, and it will therefore I think be proper to adhere to the 
1500 meters for the greatest height and depression. 

We have next to consider, what order of harmonic inequalities is most nearly 
analogous to the great terrestrial continents and oceans. The most obvious case to 
take is that of the two Americas and Africa with Europe. The average longitude of 
the Americas is between 60° and 80° W., and the average longitude of Africa is about 
25° E., hence there is a difference of longitude of about a right-angle between the two 
masses. These two great continents would be more nearly represented by an harmonic 
of the sectorial class, * rather than by a zonal harmonic, nevertheless I think the 
solution for the zonal harmonic will be adequate for the present purpose. 

Now it has been explained above that the harmonic of the fourth order represents 
an equatorial continent and a pair of polar continents. In the case of the 4th 
harmonic therefore there is a right angle of a great circle between contiguous con- 
tinents. We may conclude from this that the large terrestrial inequalities are about 
equivalent to the harmonic of the fourth order. 

Table V. (b), § 7, gives the maximum stress- differences under the centre of the equa- 
torial elevation of the several zonal harmonics, the height of each being 1500 meters* 

* The sectorial harmonic of the fourth order sin 4 6 cos 40 would well represent these two great con* 
tinents. It would fairly represent China and Australia; but would annihilate the Himalayan plateau, 
and place another great continent in mid-Pacific. It is not at all difficult to find the stress-difference 
under the centre of a sectorial inequality, but to find it generally involves the solution of a cubic 
equation. 
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The point at which this maximum is reached is given in each case, and Plate 20, fig. 4, 
illustrates graphically the law of variation of stress-difference. 

The second harmonic cannot be said to represent a continent, and the table shows 
that in each of the other cases the maximum stress-difference is very nearly 4 tons 
per square inch. The depths of the maximum point are of course very different in 
each case. 

We have concluded above that Africa and America are about equivalent to an 
harmonic of the fourth order, hence it may be concluded that the stress-difference 
under those continents is at a maximum at more than 1100 miles from the earth's 
surface, and there amounts to about 4 tons per square inch. A comparison with 
Table VII. shows that marble would break under this stress, but that strong granite 
would stand. 

The case of the isolated continent investigated in § 8 appeared likely to prove the 
most interesting one, for the purpose of application to the case of the earth. But 
unfortunately I have found it difficult to arrive at a satisfactory conclusion as to the 
proper height to attribute to the continent. 

The average height of the American continent is about 1100 feet above the sea, and 
the average depth of the Pacific Ocean about 15,000 feet. If the water of the Pacific 
be congealed into rock, it will have an effective depth of 10,000 feet. The greatest 
height of the American continent above the bed of the dried Pacific when smoothed 
down must be fully 12,000 feet or 3700 meters. The height of the great central 
Asian plateau above the average bed of the southern ocean (after drying) must be 
considerably more than this. 

Now in the application to the homogeneous planet the heights are to be halved to 
allow for the smaller density of surface rock. 

I therefore take 2000 meters as the height of the top of the equatorial table-land 
above the remaining approximately spherical portion of the sphere. 

The investigation of § 8 then shows that the equatorial table-land will give rise to a 
stress-difference of 4*1 tons per square inch at a depth of 660 miles ; and that the 
equatorial table-land counterbalanced by the pair of polar continents (the second 
harmonic constituent being absent) gives a stress difference of about 3 # 8 tons per 
square inch at a depth of 590 miles. 

This estimate of stress-difference agrees in amount, with singular exactness, with 
that just found from the case of the 4th zonal harmonic, but the maximum is reached 
400 or 500 miles nearer to the earth's surface. 

I think there can be no doubt but that there are terrestrial inequalities of much 
greater breadth than that of my isolated continent ; thus this investigation for the 
isolated continent will give a position for the maximum stress-difference too near the 
surface to correspond with the largest continents. On the other hand, I do not feel 
at all sure that I have not considerably underestimated the height of such a compara- 
tively narrow plateau. 
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In the present paper it has been impossible to take any notice of the stresses pro- 
duced by the most fundamental inequality on the earth's surface^ because it depends 
essentially on heterogeneity of density. 

It is well known that the earth may be divided into two hemispheres, one of which 
consists almost entirely of land, and the other of sea. If the south of England be 
taken as the pole of a hemisphere, it will be found that almost the whole of the land, 
excepting Australia, lies in that hemisphere, whilst the antipodal hemisphere consists 
almost entirely of sea. This proves that the centre of gravity of the earth's mass is 
more remote from England, than the centre of figure of the solid globe. 

A deformation of this kind is expressed by a surface harmonic of the first order, for 
such an harmonic is equivalent to a small displacement of the sphere as a whole, with- 
out true deformation. Now if we consider the surface forces produced by such a 
deformation in a homogeneous sphere, we find of course that there is an unbalanced 
resultant force acting on the whole sphere in the direction diametrically opposed to 
that of the equivalent displacement of the whole sphere. 

The fact that in the homogeneous sphere such an unbalanced force exists shows that 
in this case the problem is meaningless ; it is in fact merely equivalent to a mischoice 
in the origin for the coordinates. But in the case of the earth such an inequality does 
exist, and the force referred to must of course be counterbalanced somehow. The 
balance can only be maintained by inequalities of density, which are necessarily 
unknown. The problem therefore apparently eludes mathematical treatment. 

It is certain that so wide-spreading an inequality, even if not great in amount, must 
produce great stress within the globe. And just as the 2nd harmonic produces a 
more even distribution of stress than the 4th, so it is likely that the first would 
produce a more even distribution than the 2nd. 

It is difficult to avoid the conclusion that the whole of the solid portion of the earth 
is in a sensible state of stress. 

I would not however lay very much emphasis on this point, because we are in such 
complete ignorance as to the manner in which the equilibrium of the solid part of the 
earth is maintained. 

From this discussion it appears that if the earth be solid throughout, then at a 
thousand miles from the surface the material must be as strong as granite. If it be 
fluid or gaseous inside, and the crust a thousand miles thick that crust must be 
stronger than granite, and if only two or three hundred miles in thickness much 
stronger than granite. This conclusion is obviously strongly confirmatory of Sir 
William Thomson's view that the earth is solid throughout. 
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